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FOREWORD 

The  theoretical  development  of  the  determination  of  the  steady- state 
cable  configurations,  the  treatment  of  the  vortex  shedding  loading,  and 
the  work  presented  in  Appendix  B  are  due  entirely  to  Dr.  S.  Fersht. 
Programs  were  developed  by  both  Mr.  H.  L.  Butler  and  Mrs.  Mariann 
Moore.  The  program  manager  was  Dr.  A.  M.  Soldate.  Valuable 
guidance  was  given  to  the  program  by  Dr.  C.  Dudley  Fitz. 

This  program  was  accomplished  under  the  sponsorship  of  the  Advanced 
Research  Projects  Agency  (ARPA)  and  was  a  portion  of  ARPA's  effort 
to  investigate  several  problems  related  to  the  design  characteristics 
and  operational  features  of  a  high  altitude  tethered  balloon  system. 


SYNOPSIS 


A  computer  study  has  been  made  of  nonsteady  aerodynamic  loadings  on 
a  long,  cylindrical  cable  of  the  continuous  glass  fiber-resin  type  used  as 
a  tether  for  a  balloon  at  altitudes  of  approximately  100,000  feet.  No  im¬ 
portant  interactions  between  torsional,  longitudinal,  and  lateral  modes 
were  found.  Furthermore,  the  effects  of  lateral  loadings  from  gusts  or 
vortex  sheddings  were  found  to  be  unimportant.  Computer  programs 
are  presented  that  enable  computations  to  be  made  of  cable  motions 
resulting  from  localized  gust  loadings  and  from  vortex  shedding 
phenomena. 

Certain  laboratory  and  field  tests  are  recommended  for  further  studies 
of  the  effectiveness  of  the  continuous  glass  fiber- resin  cable  as  a  bal¬ 
loon  tether. 
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1.  INTRODUCTION 


Preliminary  system  analyses  have  indicated  the  technical  possibility  of 
establishing  and  maintaining  tethered  balloon  systems  at  altitudes  in 
the  order  of  100,000  feet.  The  practical  feasibility  of  such  systems  is 
critically  dependent,  however,  upon  the  characteristics  and  properties 
of  the  tether  cable.  Cable  width  must  be  kept  small  to  minimize  the 
aerodynamic  drag.  At  the  same  time,  the  cable  cross  section  must  be 
sufficient  to  carry  the  weight  of  the  portion  of  the  cable  below  or  to 
withstand  the  accumulated  drag  on  the  cable  above.  A  simple  calculation 
shows  that  constant  diameter  cables  constructed  of  ordinary  steel  are 
not  capable  of  supporting  their  own  weight  over  the  height  of  interest. 
Thus,  strength,  weight,  and  drag  properties  are  intermeshed,  and 
trade-offs  between  these  properties  must  be  made. 

If  cable  weight  were  the  only  problem,  a  possible  solution  could  be 
obtained  by  tapering  or  stepping  down  the  cable  diameter  at  lower 
altitudes  to  decrease  the  weight  to  be  supported.  Drag  force  must 
also  be  considered  however,  and  these  forces  can  only  be  balanced  by 
a  restraining  horizontal  force  at  the  ground  terminal  end  of  the  cable. 
Since  the  drag  forces  are  cumulative  (increasing  toward  the  bottom),  a 
tapered  cable  must  increase  its  diameter  at  the  lower  end.  The 
resulting  hourglass  shape  would  be  correct  for  only  one  wind  profile. 

The  hourglass  shape  is  also  difficult  and  expensive  to  fabricate  and 
to  operate. 

Similar  results  are  encountered  in  attempting  to  use  balloons  or  kite 
devices  at  intermediate  altitudes  to  help  support  the  cable  weight. 
Although  such  devices  reduce  the  load  requirements  on  the  upper  cable, 
the  drag  forces  are  increased  considerably,  requiring  an  increase  in 
cable  size.  Furthermore,  the  varying  lift  provided  by  aerodynamic 
devices  will  create  rapidly  changing  geometric  conditions,  possible 
instabilities,  and  certainly  will  increase  the  complexity  of  the  winch- 
c  able  subsystem. 
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Another  approach  for  solving  the  cable  problem  lies  in  the  improvement 
of  the  cable  materials.  In  searching  for  a  suitable  material,  a  high 
strength/weight  characteristic  is  one  of  the  first  properties  to  consider. 
Fiberglass  with  a  strength  to  weight  ratio  many  times  that  of  steel  is  a 
strong  candidate  for  the  high  altitude  balloon  cable.  Owens  Corning  has 

developed  a  combination  of  glasses  vhich,  drawn  in  individual  filaments, 

/ 

displays  a  strength  in  the  order  of  10  psi.  The  Owens  Corning 
scientists  have  also  demonstrated  their  ability  to  produce  long  multi¬ 
strand  fiberglass  cables.  Under  ARPA/NOTS  funding  in  the  fall  of 
1964,  Owens  Corning  drew  a  cable  consisting  of  30,  000  individual  fibres 
with  an  overall  diameter  in  the  order  of  0.  1  inch  and  a  total  length  in 
excess  of  80,  000  feet.  The  group  of  fibers,  bonded  together  with  epoxy, 
has  a  strength  over  the  total  area  including  the  voids  in  the  order  of 
0.  25  x  10^  psi.  Strength  to  weight  ratio  of  this  glass  cable  exceeds  that 
of  ordinary  steel  cables  by  a  factor  of  twenty.  Furthermore,  cost  is 
estimated  as  only  10  cents  per  foot.  This  constant  diameter  cable 
should  easily  carry  the  required  weight  of  the  cable  at  the  top  end  and 
resist  the  static  drag  load  at  the  ground  terminal  point. 
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Fiberglass  has,  however,  a  serious  disadvantage  in  its  weakness  to 
compressive  and  shear  loadings.  This  weakness  was  clearly  demon¬ 
strated  in  a  test  flight  accomplished  in  December  1964.  In  this  test  a 
sudden  failure  of  the  cable  occurred  while  the  cable  was  paying  out  with 
a  rapidly  rising  balloon.  Location  and  cause  of  the  initial  failure  were 
not  readily  apparent.  It  was  interesting  to  note,  however,  that  me 
sudden  release  of  tension  apparently  allowed  a  compressive  wave  to  be 
generated  which,  as  it  propagated  along  the  cable,  caused  the  cable  to 
puff  out  the  individual  filaments  at  many  locations.  The  cable  broke 
into  many  sections  and  fell  as  shards. 
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This  experience  illustrates  the  basic  need  for  an  analysis  of  the  dynamic 
characteristics  of  the  cable  and  determination  of  the  possible  occur¬ 
rence  and  magnitude  of  compressive  conditions  in  a  balloon  tethering 
cable.  From  this  analysis,  it  will  be  possible  to  gain  a  better  under¬ 
standing  of  the  required  cable  material  properties  and  dynamic  charac¬ 
teristics. 

A  complete  analysis  of  the  cable  conditions  can  be  accomplished  by  the 
following  steps. 

a)  An  analysis  of  cable  statics  and  dynamic  characteristics  at 
the  fully  deployed  condition  {balloon  at  its  float  position) 
when  the  balloon  and  cable  are  subjected  to  an  arbitrary 
wind  profile. 

b)  An  analysis  of  the  balloon,  cable,  and  cable  control  during 
launch  and  ascent  for  various  wind  conditions. 

c)  An  analysis  of  the  cable  and  cable  control  during  recovery 
operations. 

The  study  accomplished  under  the  current  program  addressed  the  first 
of  these  steps. 

As  a  result  of  this  study,  it  was  determined  that,  at  a  fully  deployed 
state,  the  cable  can  be  expected  to  be  comparatively  stable  and  that 
the  naturally  induced  vibrations  are  not  expected  to  seriously  affect  a 
cable  constructed  of  fiberglass. 

Field  tests  also  have  been  suggested  which  will  enable  substantiation  of 
the  conclusions  drawn  from  the  theoretical  investigations. 
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2.  METHOD  OF  ATTACK 

The  general  attack  used  on  the  problem  of  cable  dynamics  is  to  consider 
perturbed  motions  from  the  equilibrium  (steady  state)  shape  of  the 
tether  cable,  the  equilibrium  state  being  determined  by  the  balloon  lift, 
balloon  drag,  balloon  altitude  and  by  the  drag  along  the  cable  as  deter¬ 
mined  by  the  wind  velocity  profile. 

The  dynamic  equations  employed  are  similar  to  those  used  bv  NESCO 
for  an  extensive  and  detailed  numerical  analysis  of  the  riser  and  drill 
string  system  of  Project  Mohole  (Ref.  1).  Accordingly,  NESCO1  s 
previous  experience  in  appropriate  numerical  analysis  techniques  is 
directly  applicable. 

The  determination  of  the  steady  state  profile  will  be  discussed  first. 


3.  STEADY-STATE  CABLE  PROFILE 


Practical  considerations  dictated  the  treatment  of  cable  dynamics  in 
two-dimensional  space  coordinates  rather  than  development  of  a  three- 
dimensional  steady- state  program.  (The  treatment  of  the  steady'- state 
profile  can,  however,  be  easily  generalized  to  three  dimensions.  ) 

Geometrical  parameters  and  coordinates  are  defined  by  Fig.  1.  It  may 
be  noted  (in  this  figure)  that  altitude  is  taken  as  the  x-coordinate,  and 
the  lateral  displacement  is  indicated  by  the  y- coordinate ;  ds  represents 
an  element  of  the  original  length,  and  dS  represents  an  element  of  the 
final  length  of  the  system  such  that 

m  ds  =  mdS  (1) 

o 

where 

mQ  is  initial  mass  per  unit  length  (kg/m) 

and 

m 

m  -  is  mass  per  unit  length  (in  stretched  condition) 

with 

T 

C  1  +  j— ^  (cable  stretch  factor) 

T  =  tension  at  element 

E  -  Young's  modulus 

A  =  cross  sectional  area  of  cable 
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Definition  ot'  Axes 


PA  -  3-10.130 


Furthermore,  the  angle  between  the  element  axis  and  the  vertical  is 
designated  as  9  where 


0u 

=  TF 


Q  =  cos  6 


dx 

dS 

(2a) 

-  ll 

dS 

(2b) 

Assuming  a  horizontal  wind  in  the  (x,  y)-plane,  the  wind  loading 
normal  to  the  element  axis  is  designated  as  and 


P  =  p  Q 
rx  rn 


(3a) 


P  =  -  p 

P 

(3b) 

equations  of  equilibrium 

are  the  following. 

dx 

ds 

=  P-C 

(4a) 

ds 

-  Q-C 

(4b) 

dP 

ds 

Q-C 

T 

O’*  * 

-  PyP 

(4c) 

dQ 

ds 

P-C 

T 

P  P  - 

y 

Px  +  mg  Q 

(4d) 

-li¬ 

ds 

=  C 

V* 

mgj  P  +  pyQ 

(4e) 

Assuming  a  horizontal  wind  of  velocity  v  in  the  (x,  y)-plane,  the 
normal  wind  loading  can  he  expressed  as 

P„  =  T^cd(r>p2dvIvI  -  <5> 

where  the  Reynolds  Number  for  the  cable  is  taken  as 


R  = 

|vD  r 

1  V 

(6) 

Px  = 

p  Q 

rn 

(7a) 

py  = 

P  P 
n 

(7b) 

The  following  are 

input  parameters: 

D  =  diameter  of  cable,  meters 

EA  =  force,  kgf  -  kilogram  weight 

m  -  initial  mass/ unit  length  of  cable,  kg/m 

g  acceleration  of  gravity  at  sea  level,  9.81  m/sec^  -  used 

as  a  conversion  factor  to  express  all  forces  in  kgf 

o  mass  density  of  air  kg/m^ 

,  .  . 2, 

V  kinematic  viscosity  ot  air,  m  /sec 

Boundary  conditions  are  obtained  from  the  followin’?  v  '  i- i  •  rat  ions : 


a)  For  estimation  purposes  the  balloon  may  1:  c  considered  an 

a  sphere  filled  with  helium.  Assuming  i  cual  p  essuie 
inside  end  outside  the  balloon,  the  static  l.  i  s 


(g  is  the  value  of  the  gravitational  eoir-t  in*  at 


b)  If  the  dead  load,  which  includes  the  balloon's  membrane 
and  instrumentation,  is  W,  the  net  static  lift  is 

Wl  =  L  -  W  (9) 

where  W  is  an  input  parameter  in  Kgf;  is  an  input 
parameter  (diameter  of  balloon}  in  meters  and  1  is  the 
initial  length  of  the  cable  in  meters. 

c)  The  wind  load  on  the  balloon  is 


P(x.)  ? _ 

-  "/8 g  DbCD(R'.)v(Xi) 


v(*t> 


(10) 


in  which  the  Reynolds  Number  for  the  balloon  is  designated 


R 


v(x.) 

Db 


Cl) 


d) 


where  i/(x)  is  the  kinematic  viscosity  as  a  function 
of  altitude. 

Boundary  conditions  are: 

')  at  s  =  *  :  An  initial  guess  for  x  and  y  is  made,  then 


<d)  at  s  =  0  :  x  -  0,  y  =  Y 
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The  foliowing  data  are  required  as  input  to  the  main  steady-state 
program. 


A  table  of  wind  velocity  as  a  function  of  altitude;  the  tabu¬ 
lar  velocities  being  denoted  v(x).  The  profile  can  be  of 
arbitrary  shape.  The  particular  profile  employed  in  the 
present  study  is  shown  in  Fig.  2.  The  interpolated  values 
of  wind  velocity  are  denoted  v(x). 

Subprograms  to  compute  the  quantities  p(x)- -atmospheric 
density,  t/(x)  - -kinematic  viscosity  of  atmosphere, 
C£j(R)--drag  coefficient  of  cable  (R  is  the  Reynolds 
number),  C^(  R^  )- -drag  coefficient  of  balloon.  Atmospheric 
properties  are  essentially  those  of  the  ARDC  Standard 
Atmosphere  (Ref.  2),  the  drag  coefficients  for  a  cylinder 
and  a  sphere  are  derived  from  the  literature  (Ref.  3).  The 
expressions  for  these  quantities  used  in  the  subprograms 


Density  p(x)  in  kg/m 

pW  „  I0n/4)  -(*/ 16000) 

2  - 1 

Kinematic  viscosity  p(x)  in  rn  sec 

p(x)  is  defined  over  three  ranges  of  x 

innnn  .  .  , n-4. 823  + (x/26800) 

••  10000  meters  p  =  10 


17000  meters  v  -  10 


-5.  035  t-  (x/  17100) 


x  >  1 7000  meters 


.-5.  2421  +(x/  14250) 


Reynolds  Number  R(x)  for  the  cable 
I  V  (x)DP(s)  I 

rw  \  _ y 


ALTITUDE 


Drag  coefficient  of  cable  Cq(R) 

Cp(R)  is  defined  over  five  ranges  of  Reynolds 
number  R. 


R 

5 

2.  23 

CD 

=  10.  8  R-' 

742 

526 

R 

o 

00 

CD 

=  9.  15  R 

232 

R 

£ 

1000. 0 

CD 

=  4.95  R~* 

R 

S 

10000. 0 

CD 

=  1.0 

R 

> 

10000. 0 

CD 

=  1.5 

Figure  3  shows  the  range  of  drag  coefficients  for  the 
cable  corresponding  to  the  wind  profile  of  Fig.  2. 

For  a  cable  diameter  of  approximately  0.01  feet,  the 
Reynolds  Numbers  are  below  (by  two  orders  of  mag- 
natude)  those  for  the  transition  region  (R  <  3  x  105). 
On  the  other  hand,  the  diameters  of  balloons  are  of 
such  a  magnitude  that  transition  will  be  encountered 
at  certain  altitudes.  For  a  5-foot  diameter  sphere, 
the  transition  region  is  encountered  at  about 
50,  000  feet. 

5)  Drag  coefficient  of  balloon  Cj)(R^) 

C^(Rf)  is  defined  over  six  ranges  of  Reynolds 
Number  Rf. 


V 

1.0 

27.  4 

R~* 

961 

V 

10.  0 

s  = 

27.  4 

R'- 

804 

572 

R^ 

100.  0 

S  s 

16.  0 

R'- 

R, 5 

1 130.  0 

CD  = 

7.  25 

r'1 

4 

R, 5 

10000. 0 

CD  = 

0.  4 

Rr  > 

10000. 0 

S  = 

0.  44 

ALTITUDE  ,  I0J  ft 


Figure  3 

Drag  coefficient  as  a  function  of  altitude  for  cylindrical  cables 
and  small  spheres  corresponding  to  the  wind  profile  of  Fig.  2 
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The  method  of  solution  is  based  on  an  iterative  procedure  employing  a 
fourth  order  Runge-Kutta  technique  to  integrate  the  system  Eqs.  4a 
through  4e  along  the  length  of  the  cable.  Assuming  initial  values  for  the 
altitude  x  and  the  lateral  displacement  y  at  s  =  i,  the  tension  T 
and  the  dependent  variables  P  and  Q  may  be  determined,  thereby, 
giving  a  complete  set  of  initial  conditions  at  the  upper  end  of  the  cable. 
Values  for  the  altitude  at  the  lower  end  of  cable  are  obtained  for  various 
values  of  x  at  the  upper  end  by  the  numerical  integration  scheme.  The 
objective  of  this  procedure  is  to  choose  the  values  of  x  at  the  upper 
end  in  such  a  manner  as  to  force  the  sequence  of  values  of  x  at  the 
lower  end  xq  to  converge  to  zero  (ground  level). 

The  values  xj  and  y*  are  input  parameters  and  the  result  of  the  first 

^  i  ' 

integration  is  x  (as  well  as  y  ,  T  ,  P  ,  and  Q  ).  Succeeding  values 
.  o'  ’  o  o  o  o'  & 

of  Xj  are  determined  by 


2. 

for  i  -  1 ,  2,  .  .  .  . 

The  iterative  procedure  is  carried  out  until  the  condition 

x1  +  1  s  e  (14) 

o 

is  met,  where  e  is  an  input  tolerance. 
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The  equations  of  equilibrium  given  in  the  final  coordinate 


system  are 

dx  _  p 
dS  ~ 

(15a) 

&  -  « 

(15b) 

ai  =  t  [(p,  *  m§)Q  -  pyp 

( 1  5c) 

tfT  T  [pyP  '  (Px  *  '"*) Q  ] 

( 15d) 

df  =  (P,  +  m*) P  *  PyQ 

(15c) 

where  m  is  taken  as  m  and 

o 

S-  *  [  ('  +  EA-)ds 

(16) 

x,  y,  T,  P,  Q  at  S  -  S;  are  taken  from  the  solution  of  the  above 
iterative  procedure. 

With  the  use  of  the  wind  velocity  profile  shown  in  Fig.  2,  a  determina¬ 
tion  of  the  steady-state  cable  profile  was  made.  The  cable  properties 
were  the  following. 

Diameter  0.  002^  meters 

EA  (Product  of  Youngs  Modulus 

4 

x  cross  sectional  area)  3.451  x  10  kgf 


\ 
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E  107  psi" 

m  8. 0  x  I0*4  Kg/m 

o 

Unstretched  length  of  cable  36,  000  meters 

The  balloon  properties  were  taken  as 

Diameter  40  meters 

Dead  load  50  Kgf 

(The  properties  of  the  balloon  are  not  regarded  as  essential  to  the 
problem,  but  must  be  such  that  the  system  is  stable.  ) 

The  results  of  the  numerical  example  are  the  following:  After  five 
iterations  the  program  was  terminated  giving  the  results 


S 

X 

y 

T 

P 

Q 

m 

m 

m 

kgf 

3.  6 

x  104 

2.  75  x  104  2. 

0  x  104 

957.  1 

0.973 

•0. 230 

0.  0 

13.5  4. 

12  x  104 

746.  6 

0.427 

■0.  904 

3.  687  x 

.  io4 

The 

final  co 

mfiguretion  was  the 

n  determined  and  is 

plotted  in  Fi 

g-  4. 

Ag  ft 

•ement 

with  the  boundary  c 

onditions  for 

x  and  y 

at  the  ground 

level 

l  wa  s  w 

ithin  80  meters. 

A  fh 

nv char 

c,  listing  and  other 

details  of  th < 

e  prog r. 

un,  is  given 

in 

Appt 

•ndix  A. 

The  output  of  this 

program  is 

used  as 

input  to  the 

dynamics  program. 

^Private  conv  rsation  with  Mr.  Sheldon  D.  Elliot,  Jr.,  gives  a  value 
of  7  x  10^  psi  per  E  for  the  glass  fiber-resin  cable.  (Corresponds  to 
a  specific  gravity  of  1.  6,  ) 


* 
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STATIC  RESULTS 


HORIZONTAL  C'SPLACEMENT,  -y 


Figure  4 

Static  configuration  of  cable  with  representative  air 
veloc ity  profile 
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4.  CABLE  DYNAMICS 

Hie  primary  purpose  of  the  study  of  cable  dynamics  is  to  consider 
possible  situations  that  could  lead  to  the  existence  of  momentary 
stresses  far  in  excess  of  those  produced  in  the  static  or  steady  stale 
loading  condition.  One  possible  situation  is  the  presence  of  clear  air 
turbulence  in  a  certain  altitude  zone.  Another  situation  is  excitation 
by  vortex  shedding.  A  general  question  concerns  possible  interactions 
between  various  modes  of  cable  vibration  leading  to  stress  amplifica¬ 
tions.  The  possibility  of  the  existence  of  important  vibrational  inter¬ 
actions  will  be  discussed  first. 

With  the  coordinate  system  of  Fig.  1  generalized  to  three  dimensions 
by  a  z  coordinate  pointing  up  from  the  plane  of  the  figure,  and  the 
assumption  that  cable  is  reasonably  vertical,  it  can  be  shown  that  the 
coupling  between  torsional  and  lateral  vibration  modes  is  expressed 
by  the  following  equations  (Appendix  B). 


18 


f* 


Where  G  is  the  shear  modulus;  21  is  the  polar  moment  of  inertia;  J  is 
the  rotational  inertial  of  the  cable;  0  is  the  angle  of  twist  (torsion); 
m^  is  a  distributed  torsional  moment  per  unit  length  of  cable 
(nonexistent  in  the  present  problem),  and  p^  and  p^  are  distributed 
loads  in  the  y  and  z  directions,  respectively,  per  unit  length.  Lateral 
motions  are  coupled  to  torsional  motions  through  terms  of  the  type 

sHUS) ,nEq-  ”• 

The  importance  of  the  coupling  can  be  considered  by  rewriting  the 
7.hs  of  Eq.  17  as 


*  2L 


(?) 


*  hs  Eq.  1  7. 


00  9y  d^y  d0  dz  z 

^x  dx  ^  2  dx  dx  ~  2 
ox  ox 


2  2  2  2 

The  quantities  5  y/dx  ,  d^z/  dx  are  essentially  proportional  to  the 
inverse  radii  of  curvature  in  the  yx-  and  zx-planes.  From  an  inspec¬ 
tion  of  Fig.  4,  the  minimum  radius  of  curvature  might  be  expected  to 

4 

be  on  the  order  of  lx  !0  feet.  With  the  assumption  of  a  maximum 
shearing  strength  of  200,000  psi.  the  maximum  value  of  the  quantity 
(  *6/ dx)  is  given  by 


(assuming  G.  t  x  10  psi) 


I’he  torsional  wave  velocity  in  the  cable  is  (2G1/J)  %  1  x  1 

If  a  purely  sinusoidal  torsional  wave  is  considered,  the  ratio 
(d0/  3x)/(hZe/  txZ)  -  (A/2r),  a  being  the  wave  length. 


,  r  3  ,  -  1 

x  10  ft-sec 


-mww 


to  be  m  the  order  of  unity,  therefore,  the  period  of  excitation  must  be 
about  100  seconds.  For  such  an  excitation  period,  however,  the  situ¬ 
ation  is  becoming  essentially  that,  of  static  loading.  Unless  lateral 
excitations  can  decrease  the  radius  of  curvature  greatly  from 
the  static  values,  there  would  seem  to  be  no  effect  of  this  type  of 
excitation  upon  torsional  vibrations. 

The  effect  of  torsional  vibrations  on  lateral  vibrations  is  considered 
by  inspection  of  the  magnitude  of  the  term  2GI  (^0/dx)  which  occur' 
in  Eqs.  18  and  19.  This  term  is  a  torsional  reaction  moment  and,  with 
the  assumed  cable  properties,  is  on  the  order  of  3  foot-pounds  at  maxi¬ 
mum  shear  stress.  When  the  cableloading  N  is  considered  to  be  on 
the  order  of  1000  pounds,  it  is  intuitively  obvious  that  the  torsional 
excitations  can  have  no  effect  on  lateral  displacements. 

The  discussion  given  on  torsional-lateral  dynamic  couplings  should 
not  be  interpreted  to  mean  that  cable  twist  is  not  a  possible  problem 
caused  by  rotation  of  the  balloon.  We  assume,  however,  that  the 
attachment  of  the  balloon  to  the  cable  is  such  that  no  transmission  of 
torque  is  possible  from  balloon  to  cable  (or  from  cable  to  balloon). 

The  analysis  of  Appendix  B  (leading  to  Eq.  B-20)  shows,  in  addition, 
that  out-of-plane,  steady-state  cable  configurations,  i.  e.  ,  static 
three-dimensional  cable  configurations  will  not  contribute  to  torsion 
tf  no  end  torque  exists  at  the  point  of  cable  attachment. 
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equations  for  two-dimensional  cable  configurations  are  (using  the 
coordinate  system  of  Fig.  1) 


dS 


+  P  (t)  -  mg 

X 


(22) 


as 


+  p  (t) 

y 


(23) 


where  S  =  S(t)  is  the  deformed  arc  length  at  time  t  .  The  deformed 

arc  length  S(t)  is  related  to  the  arc  length  S  of  the  steadv-state 

o 

deformed  cable  by  the  differential  expression 


where  V  is  the  tangential  component  of  the  displacement  vector  (vector 
representing  displacement  of  a  given  point  of  the  steady-state  shape  to 
its  new  position  at  time  t  by  the  influence  of  the  dynamic  forces). 

Now 


An  -  &L 
as  AE 


(25) 
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If  AT  Z  x  10  lb,  E  -  10  psi,  and  cable  diameter  =  0.01  ft, 

(atj/a SC))  =  0.018.  Therefore,  ever,  for  a  very  large  increase  in  cable 
tension,  the  approximation  3S  %  dSQ  is  valid. 
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The  treatment  of  the  dynamic  equations.  Eqs.  22  and  23  by  the  assump¬ 
tion  of  small  departures  from  the  steady-state  configuration  proceeds 
us  follows. 


X  -  X  f  U 

o  c 


y  =  +  v 


where  uc,  v^  are  the  coordinates  of  the  displacement  vector  from 
the  steady-state  configuration  in  the  Cartesian  system  of  Fig.  1. 


Also  let 


=  To  +  Tt(t) 


Px  +  Px(t) 


Py(t)  +  Py(0 


With  the  use  of  the  steady-state  conditions  (which 
to  Eqs.  4), 


are  equivalent 


_3_  L  ±aV  o  a2x 

3so  v  o8so/  * ' me  =  mrr 


(28b) 
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and  the  approximation  5Sq  as  3S,  Eqs.  22  and  23  are  transformed  to 

,  -n  /ax  \  ,~t  dx  /3u  \  Su 

Tt  _d  (  _ o  \  aT  _ o  T  d  / _ c  \  9T  _ c  t 

as  las  /  3S  9S  as  las  i  as  as  px 

0X0/  o  o  o  \  o /  o  o 


Tt_a  fiIo\  .  t  T  j_fK\  t  yi.^c  .  t 

as  I  /  as  as  as  las  J  as  as  py 

0  0  O  O  o  x  O'  o  o  7 


With  the  assumption  that  the  terms  Tt  a/ aSQ{  Bxq/ aSQ),  etc.,  are 
negligible  with  respect  to  T  a/  aSQ  (au/aS0),  etc.,  Eqs.  29  and  30 
become 


SS  5S  3S  (T  as  )  '“x  m  2 

00  o  \  o  /  at 

2 

t  av  .  dv  „  a  v 

ai  yo  x  a  /„  c  \  t  _  c 

as  as  as  (T  as  )  py  m  2 

00  o\o/7  at 


Now,  if  the  Cartesian  components  of  the  displacement  vector  are 
rotated  in  such  a  way  that  the  components  are  tangential  and  normal, 
respectively,  to  the  steady  state  shape, 


uc  =  7J  cos  0  -  §  sin  0 


(33a) 


v  =  5  cos  0+7?  sin  0 

c 


(33b) 


w  h  e  r  e 


dx  dy 

n  _ o  .  a  o 

cos  0  -  ds  .  sin  0  dS 

o  o 
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Therefore, 


M-'OS  6  +  COS  e^"(T#)  -  sin  6J-  (T#-)  +Px 

o  o  o  oo 

l  eiia 

-  m  I  cos  o — ~  sin  v — *-  I 


(34a) 


■sin  o  t  cos 


os  e^r  (T  jp-  )+  sin  6 *T  +  pty 

O'  O'  o  o 

,)  52;  .  n5Zr?~| 

cos  d — +  sin  8  J 

bt  bt  J 


(34b) 


with  the  assumption  that  the  terms  5/ ^S^dx^/dS^),  etc.,  are 
negligible. 


Multiplying  Eq.  34a  by  -  sin  8  and  Eq.  34b  by  cos  8  and  adding 
the  two  equations, 


W  '  P1  sin  8  +  Pv  cos  8  = 

r»  '  n  '  ' 


o  '  o 


Multiplying  Eq.  34a  by  cos  8  and  Eq.  34b  by  sin  8  and  adding 
the  two  equations, 


s  i 1  ^  {  bn  \  t  t 

—  *  ~  (T  bS~  )  '  px  C°S  Stpy  s,ne  = 


o  o 


Ec|uation  33  is  the  dynamical  equation  for  lateral  vibrations;  Eq.  36  is 
the  dynamical  equation  for  longitudinal  vibrations.  Coupling  between 
the  two  equations  occurs  through  terms  depending  upon  T1  . 
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Stresses  in  the  cable  arise  from  the  force  T(t).  Accordingly,  con¬ 
sideration  is  given  to  ways  in  which  the  dynamic  component  T*  could 
become  appreciable.  An  obvious  situation  exists  when  the  balloon 
itself’  undergoes  an  erratic  movement  caused  by  turbulence- -  such 
moments  might  result  in  a  dramatic  increase  in  cable  tension  through  a 
longitudinal  signal;  however,  once  the  balloon  is  at  altitude  in  the 
neighborhood  of  100,  000  feet,  the  forces  exerted  on  it  should  be  uniform 
or  slowly  varying.  A  more  interesting  question  is  concerned  with  the 
effect  of  turbulence  on  the  cable  itself.  The  incidence  of  a  gust  on  a 
portion  of  the  cable  would  cause  increased  drag  and  a  resultant  lateral 
deflection.  If  the  lateral  deflection  is  large  enough,  the  cable  would 
undergo  a  local  stretching  which  would  then  be  distributed  along  the 
length  of  the  cable  by  propagation  of  longitudinal  and  lateral  modes. 

On  the  other  hand,  small  lateral  deflections  would  not  result  in  any 
appreciable  increase  in  cable  tension.  (In  the  same  way  that  the  vibra¬ 
tions  of  a  stretched  string  do  not  cause  an  appreciable  increase  in 
string  tension.  ) 


A  turbulent  gust  was  modelled  in  the  following  way.  The  maximum 
amplitude  of  the  gust  is  taken  as  100  ft/ sec,  varying  sinusoidally  in 
time  with  a  period  of  3  seconds.  The  velocity  distribution  in  space 
is  taken  as  a  gaussian  exponential  form,  the  width  of  which  at  10  per¬ 
cent  of  the  maximum  value  (at  a  given  time)  is  100  feet.  The  gust  is 
positioned  at  15,  500  meters  (51, 000  feet). 


The  drag  force  is  assumed  to  be  normal  to  the  cable.  Therefore, 


-p  sin  6 


(36a) 


t  t 

Py  -  P  cos  6  , 


where  p  is  the  normally  exerted  drag  force. 


( 36b) 
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I' 1-u-  drag  force  is  given  by  the  expression 


(37) 


Where  vq  =  VD(X)  is  the  steady  state  wind  profile;  v*  =  A  sinojt  •  f(x); 

A  -  100  ft/sec,  cl'  =  2ff/3  ;  p(x)  is  the  standard  atmosphere  density; 

Op  is  the  drag  coefficient  calculated  as  a  function  of  altitude  and  rela¬ 
tive  air  velocity,  and  D  is  the  diameter. 

The  dynamic  equation  for  the  lateral  deflection  of  the  cable  is,  therefore. 


(38) 


With  the  assumption  that  T(t)  =  T  ,  the  dynamics  program  presented 
in  Appendix  C,  was  applied  to  obtain  the  lateral  deflection  ~  as  a 
function  of  time,  starting  the  application  of  the  gust  velocity  A  sinci't 
at  time  t  =  0 .  The  deflections  at  the  excitation  midpoint  ('5,000  meters 
altitude)  for  approximately  3  cycles  are  shown  in  Fig.  5  .  Because  the 
dynamics  program  requires  the  solution  of  a  matrix  with  dimensions 
equal  to  the  number  of  mass  points  chosen  along  the  length  of  the  cable 
and  because  the  spacing  of  the  mass  points  was  taken  as  6  meters  in 
order  to  assure  accuracy,  it  was  found  that  the  computer  employed, 
the  CDC  3600,  was  not  large  enough  to  consider  the  motion  of  all  points 
on  the  cable.  This  difficulty  was  overcome  by  the  use  of  the  artificial 
end  conditions  that  the  lateral  deflections  from  the  static  configuration 
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were  exactly  zero  at  distances  well  remced  from  the  point  of  excita¬ 
tion.  The  points  of  the  curve  given  in  Fig.  5  were  obtained  with  two 
sets  of  end  conditions, 

a)  at  x  =  15,000  ±  2,500  meters 

-  =  0 

b)  at  x  =  15,000  ±  5,000  meters 

*  --  0 


Both  sets  of  end  conditions  gave  identical  results,  indicating  that  the 
motion  caused  by  a  localized  gust  at  15,000  meters  was  damped  out  at 
both  sets  of  end  points.  The  lateral  wave  velocity  at  the  15,000-meter 
altitude  point  was  approximately  1000  m/sec;  with  .he  relationship 
between  arc  length  and  altitude  given  by  the  static  configuration  of 
Fig.  4,  it  is  found  that  a  9-second  interval  is  sufficient  so  that  reflected 
waves  from  the  end  points  at  15,000  ±  2,500  meters  altitude  have 
reached  the  midpoint. 

In  order  to  demonstrate  the  results  of  the  dynamics  program  more 
dramatic  ally,  the  form  of  the  lateral  wave  for  various  points  in  time 
is  shown  in  Fig.  6  with  *be  end  conditions  ?  =  0  at  x  =  i  5,  000 
i  2,  500  meters. 

Because  gusts  are  statistical  and  do  not  recur  with  any  fixed  period, 
the  use  of  the  sinusoidal  excitation  term  is  somewhat  artificial.  The 
dynamics  program  is,  however,  perfectly  adapted  to  the  use  of  a 
transient  excitation.  In  the  example  chosen  here  the  first  3 
seconds  or,  indeed,  the  fir  t  1.5  seconds  of  behavior  could  be  taken 
as  representative  of  a  gust.  It  is  believed  that  the  gust  amplitude  and 
period  are  conservative. 
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Figure  6 

Lateral  displacement  along  cable  for  three  times  after 
initial  gust  excitation 
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The  displacements  shown  in  Fig.  5  are  of  the  order  of  1  meter  at 

4 

best.  With  steady  state  radii  of  curvature  of  the  order  of  10  meters, 
a  1 -meter  displacement  will  not  significantly  change  the  cable  tension. 
For  example,  assume  that  the  cable  length  increases  by  ((lO4  4  l)  - 
10  ^/!0  -  10  ,  as  a  fraction  of  the  original  (steady  state)  arc  length 

at  a  point  of  1 -meter  cable  displacement.  With  the  cable  properties 
employed,  the  cable  tension  (load)  would  increase  by  10"  x  3.451  x 

4 

10  %  3  Kgf,  whereas  the  steady-state  tension  is  approximately  800  Kgf. 


It  is  concluded,  therefore,  that  the  effects  of  gust  loading  on  a  tether  cable 
already  heavily  loaded  does  not  seriously  affect  cable  tension. 


5.  CABLE  LOADING  FROM  VORTEX  SHEDDING 


In  the  project  quarterly  report  (Ref.  4),  the  statement  was  made  that 
vortex  shedding  loads  were  not  considered  to  be  important.  A  further 
analysis  of  the  problem  was  made  and  is  presented  in  Appendix  D. 
Previous  conclusions  that  vortex  shedding  loads  are  not  important  were 
verified. 


6.  CONCLUSIONS  CONCERNING  POSSIBLE  CABLE 
LOADING  PROBLEMS 
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With  the  use  of  a  typical,  steady-state  wind  velocity  profile  (in  two 
dimensions),  it  has  not  been  found  that  interactions  between  torsional, 
longitudinal  and  lateral,  vibrations  caused  by  unsteady  wind  loading 
conditions  are  important.  Indeed,  unsteady  wind  conditions  along  the 
cable  have  not  been  found  to  cause  important  dynamic  loadings  by  exci¬ 
tation  of  any  one  particular  mode.  It  would  appear  that  balloon  motions 
are  much  more  important  than  are  direct  cable  loadings  in  causing 
cable  stress  changes.  It  has  been  shown  that  the  cable  is  weak  in  tor¬ 
sion,  and  certainly  no  coupling  should  be  allowed  to  exist  between  the 
cable  and  the  balloon  that  would  permit  a  rotating  balloon  to  exert 
torque  on  the  cable. 

It  is  entirely  possible  that  balloon  motions  and  the  weakness  of  the 
cable  material  in  compression  might  give  rise  to  destructive  effects 
in  the  longitudinal  mode.  This  may  occur  when  a  downward  motion  of 
the  balloon  is  suddenly  induced  during  ascent,  when  the  length  of  tether 
is  short.  The  cable,  normally  under  heavy  tension  loading,  is  momen¬ 
tarily  relieved  at  its  upper  end  and  an  unloading  (compressive)  signal 
is  propagated  along  the  cable  length.  This  compressive  signal  should 
travel  v.  ith  essentially  undiminished  amplitude  until  reflection  occurs 
at  some  point  of  cable  attachment,  e.  g.  ,  at  its  mooring  point.  Now 
reflection  of  a  compressive  stress  from  a  rigid  attachment  leads  to  a 
doubling  of  compressive  stress  and,  hence,  to  the  existence  of  a  net 
compressive  stress  in  the  cable.  Such  a  compressive  stress  (equal  in 
magnitude  to  the  original  tension  stress)  might  well  be  very  damaging. 
On  the  other  hand,  the  effective  point  of  cable  attachment  may  not  be 
rigid;  for  example,  termination  at  a  partially  wound  drum,  and  the 
reflection  of  the  compressive  signal  might  be  greatly  lessened  in 
s  e  v  e  r  i  I  y  . 
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7.  RECOMMENDED  LABORATORY  AND  FIELD  TESTING  OF  CABLE 
DYNAMIC  EFFECTS 

The  following  tests  are  recommended. 

a)  Laboratory  tests  on  the  effects  of  t.he  sudden  reduction  of 
tension  loading  of  the  glass  fiber-resin  cable  near  anchor 
points  of  differing  types.  Anchoring  configurations  that 
should  be  studied  include 

1)  a  solid,  vise  clamp 

2)  a  partially  wound  drum 

3)  a  mass  of  appreciable  inertia  solidly  clamped  to 
the  cable  between  the  anchor  and  the  free  end 

b)  High  altitude  balloon  field  tests  should  be  conducted  with 
special  emphasis  on  cable  motions.  Measurements  of 
cable  notions  could  be  made  by  suitably  mounted  acceler¬ 
ometers  at  intervals  along  the  cable  with  data  telemetered 
by  radio  transmitters. 

Although  the  load  of  the  transmitters  should  be  impercep¬ 
tible  on  the  cable,  this  load  can  be  reduced  by  supporting 
the  transmitters  on  small  balloons.  In  this  case,  drag  . 
from  the  balloons  would  also  be  imposed  on  the  cable. 

The  additional  balloons  would  also  allow  measurements  of 
wind  speed,  wind  direction,  and  atmospheric  conditions  at 
desired  intervals.  Thus,  motions  of  the  cable  in  known 
wind  fields  could  be  correlated. 

The  overall  configuration  of  the  cable  could  be  made  visi¬ 
ble  by  attaching  markers  such  as  flags  at  regular  intervals 
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along  the  cable  as  it  is  passed  out  during  ascension.  The 
string  could  be  photographed  from  a  distance  or  visual 
observations  could  be  made  by  telescopes  mounted  for 
triangulation. 


Cable  tension  at  the  balloon  and  at  the  ground  as  well  as  the 
cable  angle  at  these  two  positions  would,  of  course,  be 
im  portant 
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8. 


SUMMARY  AND  CONCLUSIONS 


A  computer  study  has  been  made  of  nonsteady  aerodynamic  loadings  on 
a  long  cable  of  the  continuous  glass  fiber-resin  type  used  as  a  tether  for 
a  balloon  at  altitudes  of  approximately  100,000  feet.  No  important 
interactions  between  torsional,  longitudinal,  and  lateral  modes  were 
found.  Furthermore,  the  effects  of  lateral  loadings  from  gusts  or 
vortex  sheddings  were  found  to  be  unimportant.  Computer  programs 
are  presented  that  enable  computations  to  be  made  of  cable  motions 
resulting  from  localized  gust  loadings  and  from  vortex  shedding 
phenomena.  Theoretical  results  obtained  to  date  indicate  that  the  high 
strength-to-weight  ratios  obtainable  with  the  continuous  glass  fiber- 
resin  cables  will  lead  to  an  effective  tether  for  high  altitude  balloons. 

The  present  dynamic  study  has  been  concerned  with  cable  behavior  in 
the  fully  extended  configuration.  It  is  believed  that  benefits  can  be 
obtained  by  extending  the  study  of  system  dynamics  to  include  the  motion 
of  the  balloon  and  cable  during  launch,  the  period  in  which  the  balloon 
rises  to  its  maximum  altitude,  and  recovery  operations.  Such  a  study 
should  also  delineate  the  effects  of  a  streamlined  cross  sectional  cable 
shape  and  other  techniques  for  reducing  aerodynamic  drag.  Such  a 
reduction  in  drag  might  have  an  effect  on  steady-state  configurations 
and  balloon  behavior  during  the  rise  period,  A  study  of  this  type  might 
also  better  define  the  way  in  which  the  cable  should  be  payed  out  during 
the  rise  period  (for  a  given  wind  profile)  and  would  thus  enable  wind 
performance  to  be  specified  more  effectively. 

Certain  laboratory  and  field  tests  are  recommended  for  further  studies 
of  the  effectiveness  of  the  continuous  glass  fiber-resin  cable  as  a 
balloon  tether. 
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APPENDIX  A 


STEADY-STATE  CABLE  PROFILE  PROGRAM 


Notation 


v 


♦ 

) 

) 


jr* 


k 


. 


L 


w 


D 

EA 

AMO 

VV 

G 

D1ST 

XL 

YL 

EP 


DSI 
PN 
PNP 
DSB 
PI  TER 

OPTION 


rAP 


diameter  of  cable  (meters) 
force  (Kgf^) 

mQ  -  initial  mass/unit  length  (Kg/m) 

weight  of  balloon  (Kgf) 

gravdty  at  sea  level  (rn/sec^) 

•(,  =  length  of  cable  (meters) 

» 

=  initial  altitude  (meters) 

i 

y i  -  initial  displacement  (meters) 

€  -  tolerance  defining  the  coi  vergence  of  the 
sequence  of  altitudes  computed  at  the  lower  end  of 
the  cable  (meters) 

£.S  =  arc  length  step  size  (meters) 

number  of  intervals  desired  in  final  integration 

printing  interval  (i.  e.  ,  print  every  PNP  step) 

diameter  of  balloon  (meters) 

maximum  number  of  iterations  performed  to 
obtain  convergence 

a  flag  such  that: 

1  -  spacial  history  of  cable  is  printed  after 
each  iteration 
0  -  omits  tills  printout 

a  flag  such  that 

1  -  cable  information  is  written  on  tape  for  use 
in  dynamics  program 
0  -  omits  writing  information  on  tape 


SK 


Required  Subroutines 


RUNGS 

DERIVE 

EROR 

RHOX 

NU 

CDR 


4th  order  R-K  integration  routine  for  1st 

order  system 

evaluates  1st  derivatives 

gives  error  code  and  aborts  program 

computes  mass  density  versus  altitude 

computes  kinematic  viscosity  versus  altitude 

computes  drag  coefficient  versus  Reynolds  number 
for  cable 


CDBAR  computes  drag  coefficient  versus  Reynolds  number 
for  balloon 

INTER  uses  linear  interpolation  to  obtain  wind  velocity 
versus  altitude  from  a  given  table 


Required  Data 


FORTRAN  Math 


Quantity 

Symbol 

Units 

Test  Case 

D 

D 

m 

0.  0025 

EA 

EA 

Kgf 

3.  451  x  10 

AMO 

m 

o 

Kg/m 

8.  0  x  10"4 

W 

W 

Kgf 

50.  0 

G 

H 

/  2 
in/  sec 

9.  8  i 

DIST 

l 

m 

36000. 0 

XL 

1 

m 

28000. 0 

YL 

V  Jt 

m 

20000. 0 

Required  Data  (continued) 


fortran 

Math 

Quant  it  \ 

Symbol 

Units 

Test  Case 

DT 

Not  used-leave 
field  blank 

EP 

c 

m 

10.  0 

DSI 

AS 

in 

6.  0 

PN 

PN 

bOOO. 0 

PNP 

10.  0 

PSB 

Db 

rn 

40.  0 

PITER 

6.  0 

OPTION 

1.0 

TAP 

1.0 

Table  of  \  vs.  v  (x) 
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Flow  Chart  tor  the  Numerical  Algorithm 


4 


c 


BALLOON 

ljMPUTES  STATIC  DlTllCTION  -  UNITS  IN  M-*G-StC 

DIMENSION  SSI  1000  |  ,XS(  1^00)  .YSl  iOOu)  .PS(  1GC0)  ♦•-E( 1CC0  J  .YO) .  YPtS)  » 
jnc.AO(9».V(l(JOltXV'lvO).TSS(iliOOa).lS(lt.OO)»XSSliOGG5)»PSS(1000S) 
uu  vKuN  N »NS. EN I » V  »XV  *  D»AMG .EA.G.USb 
COMMON  l\l  XSE- »  S3  •  XS  » YS 
j  ricAO  (60.3)  1  I  )  •  I  =  i  *9  ) 

>  FORMAT (9A8) 

IF  »mEAD{  1  J  )  4  .  5  A 

5  CALL  EXIT 

4  READ  (60.1)  D «  E A  * AMO .W .0 .D I c  T .XL.YL.DL.EP.DS1 . PN.PNP .DSo.P I  TER 
1  .OPTION. TAP 

1  FORMAT  (5E14.6) 

NNP-PNP 

real  (60,2)  N.IXVII )»V(I )*IS1*N) 

2  FORMAT  (  1  lo/ (6E12.6) ) 

WRI  TE  (6  ,3)  ( rit  AD  (  I  )  ,1  =  1.9) 

WRITE  (6  .9)  D.EA.AMO.w ,G.DI ST .XL .YL.DL.  £P .DS I , PN. OSB 
9  FORMAT. /5H  D  =  .1E15.5.7H  EA  =  .1E15.3.7H  MO  =  .1E15.5.6H  W  =  • 
11E15.5.6H  0  =  .U13.3//8H  DIST  =  .1E15.5.7H  XL  =  .1E15.3.7H  YL 
i-  .1E1S.5.7H  Dl  =  .  IE  IS • 5  I / 6H  EP  =  .1E1S.5.6H  DS 

3 1  =  .Iul3.-J.7n  Pi»  =  .1l15.5.6H  DSo  =  .lclS.S) 

LOGIC*  1 
oPk=DI ST/100. 

sl=dl 

<L  =  0 

I  rt.R AT  =  1 
J  T -P I  TER 
NNNP-NNP 

>j  00  TO  (11.13.14) .LOGIC 
11  LOO  I C- Loo  1 C+  1 

5  CALL  RnGX  (XL.RHO) 

CALL  Nt.  I  XL  .GNU) 

N  3  =  1 

CALL  INTER  (XL.V.XV.SVL) 

RL  =  A!33(SVL*OS.J/GNU) 

(  L  j  C  J  6  A  R  ( R  L  .  C  D  D  ) 

u  .  AS  1  :  7o  2  3 4*Ri(0«DSd*  ‘3-W 
; i  .  .j.l'o'/9Col*'.,(0*COo*oSo**2*SVL*AuS(SVL)/G 
'  ■  >.  <T  ( ‘iL * *2  +C VI. » * 2  > 

T  L 

=,iL/TL 
Y  (  1  )  =  X  L 
s'  I  2  )  =  Y  L 
r  (  3 ) =  TL 
T  1  '.  )  -  P  L 


’’  i  J  )  '=  0  L 
,  ,1  0. 


i  ■ 

i 

01  ST 

, i 

---os  I 

10 

=  \J 

V.  '• 

1.  L  ' 

'\GS 

7  '.o 

-•  I  I  ) 

-  Y  (  3  1 

■  ..  *L» 

PR 

J  - 

1 

‘  ... 

(  J  )  - 

3 

>:  . 

(  J  )  - 

Yin 

y  j 

(  J  )  - 

Y  (  2  ) 

r  ■> 

(  J  )  = 

Y  1  3  ) 

1  J  )  - 

Y  ( ■'»  ) 

( S.US.S.Y.YP.ID) 


1 
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I 


o 


V- 


uS' J)=  Y ( & } 

o  (.Al  L  Rl'NQS  IS.DS.5.Y.YP.ID) 

<  =  .♦1 

rss ( I  i=r i 3i 

IF  (OPTION)  91.92.91 

91  IF  ID1ST-S-GPRS)  92.95.95 
95  DPRS  =  DPRS-*0P3 

J:j*i 
SS( J 1 =  S 
X5IJ)-  Yill 
YS(J)=  Y(2) 

T  S  (  J  I  =  Y I  3 1 
P£  (  J I  =  V(4) 

0S|J)=  Y I  5 ) 

92  IF  IS-. CCCll  7.7.6 

7  IF  luPTIOM  94.50.94 
94  JJ*J*1 

SSI  JJ )  s  s 

XSI JJI  -  VI 1 ) 

YSIJJ)  s  YI2) 
rsu-'l  =  Y  I  3  I 
PSIJJ)  =  Y I  4 ) 

CSI  )  -  Y  15  I 

h'IT.  (6  .IB)  laSI 11  ).X5< I  I ) .YSI I  I ) »TSI 1  I ) .PSI I  1  ) .uSI I  I ) ,i I -l.JJ) 

go  re  so 

13  LCoiC  =L0G1C* 1 
Yl=AcilYll)) 

2'»  lUR.'.T-iTtRAT*l 

IF  IVI1))-£P)  20.20.86 

86  IF  (IT'RAT-JT)  26.d7.87 
26  XL=XL-Y( l)/2. 

..-0  10  S 

1 4  IF  I AuSI Y  I  1 1  l-Yl )  15.16.16 

15  Yl-At:S(Yll)) 

00  10  24 

lo  al*XL*Y1/6. 

4 1.  -  <  L  ♦  1 

T O'  I8,54),KL 
."ii  oT^Rll) 
iJ  /  ....  I  T  c  .6  .od  ) 
u':  I  .  ■>  *A  T  (  1  Iff  1  D  I  VtROLNCE  ) 

(J  "lit  (6  »  1  o )  L>lbT*XL»YL»TL*Pl»uL.S.(YIJ)*J=l*5) 

lo  1 vRVAT  I///9X.1HS.18X,1HX,18X.1HY,18X.1HT.18X,1HP,18X,1H0// 

1  I  6c 19.7)  ) 

-  1  •  ♦  I  TSSI  ll  +  TSSl  I  )  )/l  2.*EA) 

II  *  1-1 
DO  4  *  J  =  2 » I  I 

4-  SO''  -  5  *.'-*♦■*  i  •  +  TSSI J I /LA  ) 

5  >'■'  =  SUM  *  OS  1 
*  'I  Tc  (6  ,65)  SUM 

■-5  ‘'-'•'AT  (///  16H  TOTAL  LENGTH  =  .1E18.6) 

>1-  1  YL  -YI2) 

S  =  SUM 
OS  =  -SUM/PN 
■'•5  =  1 
!  =1 
J  =i 
Y I  1  ) =XL 
Y  I  2  )  —  Y  L 
Y( 3 ) =TL 
Y (4 ) =PL 
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V.' 


L 


■rJ 


f 


Y(3)=0L 
SSI1I  =s 
XSll)  =XL 
YSC 1 |  =YL 

rsil)  =TL 
Psm  =pL 
3SI1I  =QL 

ent*i. 

10-0 

XSSI  1 !*YI  1) 

TSsi 1 )-Y(3) 

PSSI  1)=Y(4) 

CALL  RUNGS  (S»DS»5»Y»YP,1D) 

60  CALL  RUNGS  I S t OS  *  5  * Y . YP . 1 D ) 

Ul  +  1 

XSS I ! 1 s Y I  1 1 
TSS  I  1  )  =Y  I  3  1 
PSSI  1  )  »  Y  (  4  ) 

IF  I  1  -NNNP 1  61.62.62 
62  NNNP-NNNP+NNP 
J  =  J  +  1 
SSI J)=S 
XS(J1=  rill 
YS I  J ) =  Y I  2 1 
TSI J  1  =  Y I  3  ) 

PS(JI=  Y ( 4 1 
v.S  I  J  1  =  YI51 

61  IF  IS*. 00011  64,64,60 
64  JJ  =  J+l 

SS(JJ)  =  S 
XSI JJ)  =  Y 1 1 ) 

YSIJJ)  =  Y I  2 1 
TSI JJ1  =  Y I  3 1 
PSIJJ1  =  YI41 
OSI JJ)  =  Y I  5 1 

WRITE  (6  .16)  tSSIL),XS(L),YSIL)»TSlL),PSIL)»uSIL),L-l,oj) 

IF  (TAP)  ’.63,1000,163 
163  WR  1  IE  i  6  •,  8  0  0  0  )  1  -DSI  .01ST 

8  0  Ju  FORMAT  I lHl , 1 9, 25H  POINTS  AT  A  INCREMENT  0F,E15.8,12H  STARTING  AT, 
1E15.8)  . 

.'.RITE  (2)  (XSSIU  . L  =  1  , 1  ) 
ftRlTt  12)  I TSSIL) ,L=1  ,1  ) 

.v RITE  12)  (PSSIU  ,L  =  1,S) 

GO  TO  1000 
END 

SUoRGUTlNE  RhOXlX.RHO) 

C  MASS  DENSITY  VS  ALTITUDE 

RHO=,25-X/16000. 

RHO=10.**RHO 

RETURN 

END 

GUbROUTlNt  NU  (X«GNU) 

C  UNEKAT1C  VISCOSITY  VS  ALTITUDE 

IF  (X-1000C. )  1,1,2 

1  GNU  =  -4.b2  J  •■X/26B00. 

GO  TO  5 

2  IF  (X-17700. 13,3,4 
b  u N U = X /  1  7K0.  -6 • 03  S 

oO  TO  6 

..  ONU  =X / 1 4250 •  -  6  •  242 1 
6  GNU=10.*«GNU 


1 

} 


i 

I 

1 
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RETURN 

END 

SUBROUTINE  INTER  I  XL  <  V .XV  < S VL ) 

OIVtN  Tn£  ALTITUDE,  THIS  ROUTlNt  INTERPOLATES  LINEaRly  TO 

ucTAIn  The  «INU  vcLOCITY  FROM  INPUT  TAoLE 

DIMENSION  VI  I0UI  .AVI  10UI 

COMMON  N#NS 

NN-N-NS+1 

IF  IXL-XV(NN))  5,5,6 
6  NS  =  1 

5  00  1  1 =  NS*N 
JsN-H-1 

IF  I XL-XV ( J ) I  1,3,4 

3  SVL  =  V ( J I 
00  TO  2 

4  SVL=(V< J+1I-VU1 )»(XL-XV( J+1))/(XV< J+lI-XVI Jl l+VI J+ll 
00  TO  2 

1  CONTINUE 

2  NS= J+ 1 

ke  turn 

END 

subroutine  cobar  <RL,cOdi 

LRmG  COEFFICIENT  VS  REYNOLDS  NUMBER  FOR  BALLOON 
IF  ISL-l.C!  1.1,2 

1  03  =  27. 4*Rl**<  .961  I 
f-0  TO  10 

2  IF  I R  L “  1 0 . I  3,3,4 

i  Cjc=27.4*Rl»»(-.&04) 
cO  TO  10 

4  IF  (RL-100.)  5,5,6 

5  C0d=16.0*RL**(-.572 I 
GO  TO  10 

6  IF  (RL-1130.)  7,7,8 

7  CDB=7.25*RL»»(-.4I 
00  TO  10 

8  IF  (RL-’OUOG.I  9,9,11 

9  CDB=.4 


uj  iu  io 
11  C  Dt3  =  •  44 
IU  RtTURN 
t  NO 

SUBROUTINE  RUNGS  (X,H»N»Y»YPRIME, INDEX) 
x  independent  variable 

H  increment  delta  x,  may  be  changed  in  value 

N  NUMBER  OF  EQUATIONS 

YPmMF^BW  ^RMO,i  bL0C^  0NE  DIMENSIONAL  ARRAY 
YPRIMt  DER I V«  T I v  t  CLOCK  uNt  DIMENSIONAL  ARRAY 

THE  PROGRAMMER  MUST  SUPPLY  INITIAL  VALUES  OF  YI1)  TO  YIN! 
INDEX  Is  A  VARIABLE  WH I CH  SHOULD  BE  SET  TO  ZERO  HFFORF  FAru 
IB.I1AL  E«m»  10  I  He  SUMOUIlHil  !  !o  MM  ,  ™«  S 
StT  OF  EQUATIONS  OR  TO  START  WITH  NEW  INITIAL  CONDITIONS. 

THl  PROGRAMMER  MUST  WRITE  A  SUBROUTINE  CALLED  OtRIvt  WHICH  CUM 

futls  the  derivatives  and  stores  them 

SUBROUTINE  DLRIVEIX.N.Y.YPRIME) 

1  DO  2  1  =  1, N 

will* =H»YPR I  ME (  I  I 

l  Z(II=Y(I)+(wl(I)*,5l 

A  =  X .n/ 2 . 

CALL  OCR  |  VL  (  A  .  N  »Z  »  YPP.  IMt ) 


DO  3  1=1. N 
til  ill  =H«VPR1ME(  1 1 
3  Z<  1  I =Y( 1 I  +  *5«W2( 1  I 
A=X*H/2» 

CALL  DERIVE! A, N.Z.YPR1ME  ) 

00  A  1=1.N 
„3( j  J=H«YPR1ME(  1  ) 
it  n  n=Yi  1 1  +w3 1  i  i 
A  =  X+H 

CALL  DERIVE  (A.N.Z.YPR1ME) 

DO  7  1=1 tN 

A'4(  1  >=H*YPR1ME(  i  i 

7  Y(1)=Y(1|+(((2«*(W2<1>+W3(l)l)+Wi(l) +W4<  1 11/6.  ) 
X=X  +  H 

CALL  DERIVE  (X.N.Y.YPR1ME) 

GO  TO  6 

5  CALL  DERIVE  ( X . N. Y  .  YPR1 ME ) 

1 NDEX= 1 

6  RETURN 
END 

SUBROUTINE  DERIVE  (X.NN.Y.YP) 

DIMENSION  VI 100) *XV( 100) *Y( 51 »YP (51 
CUMMuN  N.NS.ENT.V.XV.D.AMO.EA.G.DSB 
Z  —  Y  C  1  > 

CALL  RHOX  (Z.RHO) 

CALI.  NU  (  Z  i  GNU ) 

CALL  INTER  (Z.V.XV.SVL) 

R  =  AQS<  SVL*0*Y< 41/GNU) 

CALL  CDR  (R.CD) 

PN  =.5«RHO*CD*Y( 4)**2*D*SVL*ABS( SVL ) /G 
PX  =  PN*Y I b I 
PY  =  -PN*Y(4) 

IK  (EM)  2.2.1 
2  DEM =  1.+  Y ( 31/EA 
AM*  A'-'O/DEM 
PGM  =  PX  +AM*G 
YP! 1 ) = Y ( 4  1  *DEM 
YP ( 2  )  =  Y(5|*DEM 
YP ( 3  )  s  (PuM»Y(4)+PY*Y (SI I  *  D  E  M 
YPIM  =  YP  (-2  )  *  (  PGM*Y  (  5  )  -PY*  Y  (  4  )  )  /  Y  (  3  ) 

YP<  b  )  =  Y P (  1  |*(PY*Y (41 -PoM»Y ( b)  l/Y ( 3) 

KtTURN 

1  PGM  =  PX  .AMC*G 
YP  (1)  =  1(4) 
y  P  (  2  )  =  Y  (  b  ) 

YP (  3  I  =  PGM* Y ( 4 )  +  PY* Y ( b ) 

YPI4I  =  (  PGK»Yl  5) -PY*Y  (4  |  )*Y(  b  )/  Y  (  3  ) 

YP ( 5  )  =  ( P  Y * Y ( 4  )  - PGM*  Y ( 3 ) ) * Y ( 4 ) / Y ( 3 ) 

RETURN 

END 

SUBROUTINE  CDR ( R i CD ) 

DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER  FOR  CABLE 
IF  (R-2.23)  1»1». 

1  CD= 10.6*R»» ( -.742 ) 

GO  TO  10 

2  IF  <  R  —  8 •  0  )  3.3.4 

3  CD=9.1S»R»»(-.S26) 

GO  TO  10 

4  IF  (R-1CC0.0I  5.5.6 

5  CD=4.95*R**(-.232) 

GO  TO  10 
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I 


i 


o  If  (R-)OCOO.)  7.7.8 

7  CD=1.0 
00  TO  10 

8  C0=1.15 
lo  RETURN 

END 

i,  J5R0UT 1 Nt  EROR  (  1  I 
I  *1 

WRITE  (61.11  1 

1  FORMAT  (13H  ERROR  CODE  *.1151 
CALL  EXIT 
RETURN 
END 
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i 


I 


Y‘ 


5*'* 


| 


APPENDIX  B 

SIGNIFICANCE  OF  TORSION 
by  D  r.  S.  Fe  rsht 


I 


Consider  a  three-dimensional  right  handed  Cartesian  system  of  co¬ 
ordinates  (x.  y,  z),  fixed  at  the  top  of  the  riser  with  x  (vertical)  directed 
downwards.  In  each  horizontal  section  of  the  riser  we  have  three  forces 
and  three  moments  in  the  direction  of  the  axis  (Fig.  B-l). 

Considering  an  element  of  the  riser  there  are,  in  addition  to  the  internal 

forces  and  moments,  external  loacis  and  inertia  forces  in  the  y,  r. 

directions.  Finally,  a  distributed  external  moment  m  will  be  consid- 

x 

ered.  Coupling  of  motion  in  the  "x"  direction  with  "y",  "z'1  motion  is  a 
separate  problem  which  is  not  considered  here,  so  there  is  no  equation 
for  "x"  motion. 


Deflections  due  to  shear,  and  rotational  inertia  about  axes  normal  to 
the  deflected  shape  of  the  riser,  are  factors  which  need  be  considered 
only  when  there  is  reasot  to  expect  deformations  in  modes  whose  wave 
lengths  are  in  the  neighboi  hood  of  the  diameter  of  the  riser.  For  the 
riser,  major  exciting  forces  with  periods  lower  than  0.003 
seconds  would  be  required  before  the  inclusion  of  the  shear  and  rotary 
inertia  would  lie  warranted,  and  those  complicating  factors  are  therefore 
neglected. 


file  remaining  equations  of  motion  for  the  elemental  Fig.  B-l  are 
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Eliminating  S  and  S  from  Eqs.  B-2,  B-3,  and  B-4, 
V  z 


'M  ,  ~M 

x  y  \ 


-M 

•z  z 


~x 


X  -  X  "X  -  X 


=  J 


>  *  f?)2 .  (a 2 

V*X/  \~X/ 


M2 


(B-S) 
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Since  1  \^xJ  "r  t*le  ^ast  equation  can  be  written 
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We  now  find  the  stress  strain  relations  for  the  beam.  In  order  to 
understand  the  geometry  of  deformation  of  such  a  beam,  and  derive 
the  stress  strain  relations,  we  have  to  use  some  of  the  concepts  of  a 
space  curve.  In  our  case  the  problem  is  m  icn  simplified  by  the  fact 
that  the  cross  section  of  the  riser  has  circular  symmetry.  This 
permits  each  point  along  the  axis  of  the  riser  in  the  deformed  shape  to 
be  described  by  a  radius  vecto’- 

r  "•  >:  i  *  y  j  +  z  k  ( B-  6) 

where  l.  j  and  k,  are  unit  vectors  along  the  axes  of  the  coordinate  system. 
Assuming  that 
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the  unit  vector  tangent  to  the  deformed  riser  is 
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Hie  plane  in  which  the  riser  lies  at  each  point  is  ti.e  osculating  plant. 
This  plane  is  determined  by  two  vectors;  the  vector  t  and,  a  unit  vector 
n  which  is  normal  to  the  riser.  It  is  known  from  differential  geometry 
that 
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(B-8) 


where  p  is  the  curvature  of  the  deformed  riser,  measured  in  the 
oscillating  plane.  The  unit  vector  normal  tc  t  and  n,  i.e.,  the  unit 
vector  normal  to  the  osculating  plane  is, 


b  =  t  x  n  =  p 
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Elastic  relationships  for  bending  and  torsion  are 


P 


Mb 
El  ’ 


£0  .  T 
2GI 


(B- 101 


By  ordinary  sign  coir  ions,  the  bending  moment  vector  for  a  cross 
section  with  an  outer  normal  t  ,  is  in  the  direction  of  b,  while  the 

torsion  moment  vector  is  in  the  direction  of  t.  Using  Eqs.  B-7,  B-8, 
B-9  and  B-10,  one  can  write 


—  FI  — 

M,  b  =  —  b 
b  P 


,2  .2 

El  — Vj  +  El  — L  k  (B-  1  1  ) 

-x  -  X 


T  t  =  2  Cl  t—  t  =  T  i 
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The  sum  of  the  components  of  these  two  vectors  alonu  the  axes  are, 
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Using  Eqs.  B-l,  B-3,  B-4,  B-5,  B-10  and  B-13,  one  obtains 


(B-16) 


The  last  term  on  the  right  side  of  Eqs.  B-15  and  B-l 6  is  small  compared 
to  the  other  terms,  and  may  be  neglected. 


We  now  have  differential  equations  for  the  three  unknow  functions, 
y(x,  t),  z(x,  t)  and  0(x,  t).  For  the  static  case  this  system  of  equations 
can  be  reduced  to  the  form 
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These  equations  of  statics  may  be  used  to  gain  quantitative  knowledge 
about  coupling  between  bending  and  torsion.  Let  us  begin  with  the 
effect  of  out  of  plane  bending  on  torsion.  Assuming  that  m^  =  0,  one 
can  see  that  the  last  two  terms  in  Eq.  B-17  represents  the  bending  effect 
on  the  torsion.  Equation  B-17  can  be  readily  transformed  into  the  form, 


£[*♦©  ♦&) H &[*♦(&)  ♦(*)]■<> 

Dividing  Eq.  B-17a  by  T  j^l  +  +  (clx^  J’  °ne  °^ta^ns 


j_  dT  1 
T  dx  2 


*♦(£) 


d_ 

dx 


0 


■ ♦  (£)‘ 


(£)2]  - » 


Integrating  this  equation  gives 

T  ■  to[‘  +  (£)2  +  (s)2J 


1/2 


(B-  1 7b) 


(B-20) 


Thus,  for  the  small  angle  theory  assumed  throughout  and  verified  by 
calculations,  the  torsion  in  the  beam  does  not  vary  with  "x"  unless 
there  are  external  moments,  m^,  applied  continuously  or  discretely 


along  the  length  of  the  riser, 
torsion  loads  or  stresses. 


Out  of  plane  bending  does  not  create 
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APPENDIX  C 


CABLE  DYNAMIC  PROGRAM- -EFFECTS  OF  LATERAL  GUST  LOADING 


1 

i 


j 


Equation  of  Motion 


±(r  b*\  =  m  t. 

dS  \  /  o 


2p  p(x  ) 


-t-^DCt^(R  )p‘- 

2g  D  m  m 


-  V 

V 

o 

o 

m 

m 

— 

where 

f  2.3(x  -  x  ,j\2'| 

t  /.  A  )  \  m  midy  V 

v  =  ( A  sin  ojti  exp  < - i > 

l‘/^Al-*A2)V 

and  XA2  b  cund  the  forcing  function  v*,  otherwise  vt  =  0,  and 
Xmid  *s  tbe  cab*e  se8ment  midpoint. 


Difference  Analogue 


( ~m+ 1 ,  n+1  +  ^  ^m+1,  n  +  ^m+1 ,  n-  1  "m,  n+  1  "  m,  n 

i.n-l)  ^m-1  (^m,n+l  +  ^  ^rn,  n  +  ^m,  n-  1  ^m-l.n+l 


^ ’m-  1 ,  n  ^m-l,n-l)j 


m  f  ~ 

(5m,n+l  •  2?m,n  *  A.n-l)  ‘  DCD(Rm,Pm(  ' 

(  ' nil  1  ‘m,  n- 1  \  t  m,n  3m,n-l  \ 

A  ^  )  V°.n  '  A  S*  ) 


V  +  V 

o 

m 


v  v 
o  o 
m  m 


V  +  V 

o 
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Boundary  Conditions 


’  1 , n+ 1 


'M  •  n+1 


-  0  ,  where  m  =  1,  2,  .  .  .  ,  M 


V 


&*  r 


Initial  Conditions 


§  n  =  ° 
m )  u 


where  m  =  1 , 2 ,  .  . .  ,  ivl 


Method  of  Solution 


If  A  ,%  ,+  B  ,?  ,+C  ,  ,  £  ,  =D  for 

m- 1  m-l,n+l  m-1  m,  n+ 1  m- 1 ,  n+ 1  '  m+ 1 ,  n+ 1  m-  1 

m  =  2, ...  ,  M- 1 ,  then 


At2  Atp(x  )DC  n(R)Pl 

B  =  1+ — =i - (T  +  T  ,)  + - - - 2 — 21 — 2L 

m’1  4AS2m  m  m-2  4Smo 

o 


,  t  /  'm 
v  +  v  -  - 

°  V 

m  \ 


,  n  ’m, n-  1  \ 

~) 


D  i  -2'  -  5  ,  -  C  ,  (2*  +  §  ,  .  -  2%  -  =m, 

m-1  m,n  m,n-l  m- U  m+l,n  m+l,n-l  m,n 


,(2? 

1-  1\  m,  n 


+  r  ,  -  2?  .  -  ? 

m, n- 1  m-  1 , n  m 


1 ,  n-  l) 


Atp(x  )DC  (R  )Pe‘ 

_ m  D  m  m 

4gm  v  o 

o  m 


3  1 

m,  n-  1 


At2p(x  )DCT.(R  )P2 

r  m  D  m  m  v  v 


L.KU1 

6  O  _ 

/  '  m,  n  ^m,  n-  1  \ 

V  * - ) 

J 


V  -  (v  +  V  J  I  V  +  V 

o  o  to  / 1  o 

mm  \  m  /  I  m 


MWi  ***#  *** 


The  coefficients  A  ,  B  ,  C  and  D  define  a  set  of  M-2  linear 
m  m  m  m 

equations  for  the  unknowns  (m  =  2...,  M  -1).  These  coefficients 

are  expressed  in  terms  of  the  previous  displacements  n 

E  which  are  assumed  to  be  known.  The  matrix  of  the  linear 

m,  n-2 

system  for  the  £  takes  the  form  of  a  three  term  M-2  by  M-2 
7  m,  n  7 

band  diagonal  matrix.  This  matrix  is  triangularized  by  use  of  gaussian 
elimination  techniques  and  is  readily  inverted,  thus  yielding  the  solu¬ 


tion  for  the  £  .  The  solution  for  the  £  is  thus  carried  out 

■  m ,  n  m,n 

point-by-point  in  time,  beginning  with  the  initial  conditions  for  E^ 

i.  e.  ,  £  =0.  In  other  words  at  each  time  n^t,  the  complete  set 

m,  o  r 

n  (m  =  2.  . .  ,  M  -1)  is  obtained  and  then  together  with  previous 

5  .  etc.  ,  is  used  to  solve  for  the  complete  set  £  , 

m,  n  - 1  'm,  n  41 


:':(The  subscript  m  labels  the  space  net;  the  subscript  n  labels  the 
time  net. ) 
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Following  the  computation  of  the  above  sequences,  the  solution  for 
§(m,  t)  may  be  written  for  m  =  M  -  2 ,  M  -  3, ....  2  as 


^M-l.n+l  °M-1 


m,  n+  1 


D  -  C  ? 
m-  1  ^m- 1  m+ 1 ,  n+ 1 


Required  Data 


a)  Input  constants  .  mQ,  g,  D 

b)  Static  results.  X(S),  T(Sh  P(S)  at  intervals  As,  cable  segment 
endpoints  Xj(S)  and  x^(S)(Xj  >  x^}  and  endpoints  of  forcing 
function  range  xA1(S)  and  xA2(S)  (xA1  >xA?)* 

c)  Supporting  subprograms  to  compute.  vq>  p,  C^,  and  V 


;Note:  For  the  remainder  of  the  discussion  X  will  be  used  in  place 
of  x  for  the  altitude  coordinate. 


60 


Card  Input  to  Program  Balloon 


Columns  10 1 

20 

3C 

Card  1 

NDA 

NMO 

NYR 

Card  2 

NTAB 

NXPRT 

NVEL 

Card  3 

IDEBM 

IDEBB 

IDEBA 

Card  4 

XI 

X2 

EPS 

Card  5 

DS 

G 

D 

Card  6 

A 

XAl 

XA2 

Columns  12 

24 

36 

48 

60 

Card  Group  7 

XV(1) 

XVfNVEL 

V(  1) 

V  (NVEL) 

XV  (2) 

V(2) 

XV(3) 

Tape  Input  to  Program  Balloon 

Logical  tape  Z  contains  3  binary  records. 

a)  Record  1.  X(l) . X(NTAB) 

b)  Record  Z.  T(  1) . T(NTAB) 

c)  Record  3.  P(  1),  .  . .  ,  P(NTAB 

Nomenclature 

NDA,  NMO,  NYR  date 

NTAB  number  of  data  points  tabulated  by  static  program 

NXPRT  increment  in  X  to  use  in  printing  the  cable  section 
results  (i.e.  ,  print  every  NXPRT  x) 


- - - - -  - 


\ 


NVEL 

IDEB 

XI 

X2 

EPS 

TTEST 

DTPRT 

DTMPRT 

DT 

DS 

G 

D 

EMO 

XB  ALL 
SB  ALL 
OMEGA 
A 
XA1 
XA2 
XMID 


number  of  entries  in  wind  velocity  tables 

card--leave  blank--used  for  debugging 

X  value  at  upper  end  of  cable  segment,  meters 

X  value  at  lower  end  of  cable  segment  XI  >  X2,  meters 

€  =  tolerance  for  deflection  at  ends  of  cable  segment 

time  to  test  the  program  and  halt  at  the  first  min  and 
max  after  t  =  TTEST 

increment  in  t  to  use  in  printing  the  cable  section 
results 

increment  in  t  to  use  in  printing  the  time  history  of 
the  cable  section  midpoint 

At  =  increment  in  t  to  use  in  evaluating  the  differential 
equations 

AS  -  increment  in  S  at  which  the  static  date  was 
recorded,  meters 

.  2 

gravity  at  sea  level,  meters/sec 
diameter  of  cable,  meters 

=  initial  mass  per  unit  length  of  cable, 

Kg/meters 

x-coordinate  of  balloon,  meters 
length  of  cable,  meters 

u)  ■-  angular  frequency  of  forcing  function,  radians 
amplitude  of  forcing  function  v*,  meters/sec 
upper  limit  of  forcing  function  v*,  meters 
lower  limit  of  forcing  function  v" ,  meters 
cable  section  midpoint,  meters 
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PROGRAM  BALLOON 


J  l 


\ 


m 


PROGRAM  BALLOON 

C  DYNAMIC  ANALYSIS  uF  A  2-0  TETHLReD  dALLOON  SYSTlM  SNJ44-712 
DIMEN  S I  C'N  SPACt  I  10000)  »xl  (20u0»3  )  .TTAd!  2000  )  »PT  Ad  I XOOO ) 

1  »XT  AB ( 2000 ) 

DIMENSION  XVI 100) .VI 100) .aIMIDI 10C0 ) , STUFP ( 1 0000 ) .STUF T 1 10000 ) 
COMMON  SPACE. XI .TTAb.PTAd. I DEbM, I DEdB » I DEbA, IDtdG. IDEBS 
l.XIAB 

COMMON  /l/  STUFP, STUFT 

C  N  POINTS  DESCRIBE  INPUT  CAdLE  SECTION  otTwEEN  XI  AND  X2 
C  WHERE  XI  MORE  THAN  X2. 

C  CAdLE  ANALYZED  FROM  XI  TO  X2  AT  EACH  DT 

C  THE  CABLE  IS  SBALL  LONG  AND  HAS  NT  Ad  INPUT  INFORMATION  POINTS 

C  AT  EVERY  DS  POINT  AND  IS  RECORDED  ON  TAPE  UNIT  2 

C  NVEL  IS  NO.  OF  PTS.  IN  VELOCITY  TABLE  OF  ALTITUDE  XV  VS  VELOCITY  V 

C  IDEB  INDICES  PRINT  INTERMEDIATE  OUTPUT  IF  NOT  ECUAL  TO  ZERO 

C  DISPLACEMENT  IS  PRINTED  AT  EVERY  NXPRT  X  VALUE  IN  TIME  INCREMENTS 

C  OF  DTPRT 

C  DISPLACEMENTS  AT  THL  MID  POINT  OF  THE  CAdLE  WILL  BE  PRINTED  AT 
C  TIME  INCREMENT  DTMP.RT 

C  WHEN  T  =  T1 EST  THE  PROGRAM  IS  STOPPED  AT  THE  FIRST  MAX  AND  MIN 

C  AFTER  TTEST 

C  1ST  DISPLACEMENT  AFTER  Si  AND  LAST  BEFORE  S2  MUST  BE  WITHIN  EPS 

C  V  TO  THE  T  =  A  •  Si  N  ( OMEGA  •  T  )  '.  EXPONENTIAL  FUNCTION  OF  X 

C  SPACE! 1-2C00)=A,  l 200 1-4000 ) =  d ,  I A00 1-6C0G ) =C  .  ( 600 1 -8000 ) =D 

C  SPACE! 8CG1-1C000) =  D.RHOI XM ) .CD! KM )P ( XM) / ( 2G I 

1  READ  15.6000  NDA.NMO.NYk 
IF  ( NDA )  99.99.2 

2  READ  15.6000)  NTAU. NXPRT, NVEL 

READ  ( 5 .8000 ) IDEBM. IDEdd, I DEdA .  IDEbG . I  DEBS 
800U  FORMAT (7) 10) 

READ  (5.6001  1X1.X2  .EPS.TTEST.DTPRT.DTMPRT.DT.DS.  G.D.EMOvXdALL  . 

1 SBALL. OMEGA 

8001  FORMAT  (7F10.8: 

Read  (5.8001)  A.XA1.XA2  ,XMID 

8002  FORMA T ( 6E 12. 6  )  * 

READ  (5.8002)  ( XV ( I  )  ,V( I ) , I  =  1 .NVEL ) 

C  INPUT  X.T.P 

REWIND  2 

WRITE  ( 6 .905  0 ) NT  AS 
9050  FORMAT  I//.53X.7H  NTA6  =  » I  6  . / /  ) 

READ  ( 2) (SPACE! I i  .1 =1»NTAB> 

READ  ( 2 ) ( STUF T ( I ) .1=1 *NT AB ) 

READ  ( 2) (STUFP! I )  ,1-l.NTAB) 

DO  30  I  =  1 «N  T  Ad 
INDEX=  I 

IF  (Xl-SPACE! I )!30,AQ,4C 
3^  CONTINUE 
<tu  DO  5C  I  -  I  NOE  X  .NT  AB 
NXA1= I 

IF  I XA1 -SPACE (  I  i  )  50.60,60 
5'J  CONTINUE 
bJ  00  62  I=NXA1,NTA6 
NM I  0=  I 

IF  (XM ID- SPACE (  I  )  )  62  ,65,65 
62  CONTINUE 
65  00  70  I =NMI D »NTA3 
NXA2- I 

IF  ( XA2-SPACE  !  I  )  )  70,60,80 
CONTINUE 

oj  DO  90  !=NXA2,NTA8 
LAST= I 
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If  :  ,2  - 5. PACE ( I  )  I  90.95.95 
9 v  CONTINUE 
95  N  =  L  A  S  T - 1 NDE  X  + 1 

NXA 1=NXA1-INDEX+1 
NXA2-NXA2-I NDEX+1 
NM I D  =  NMI 0-1 NDEX+1 

WRITE  (6.70C0)  NTA5 . N. I NOEX .LAST , NX A1 ,NX A2 .X TAo ( 1 1 . TTAo ( 1 1 »PTAb( 1 ) 
1.  XTABIM .TTA6IN) .PTA6IN) 

7vJj  FORMAT (61 10./.6E20.8) 

WRITE! 6.9000 INDA.NMO.NVR 

9000  FORMAT!  lHl,30X.A5r(DTNAMlC  ANALYSIS  OF  A  T  t  THF.RED  BALL-SN344-712. 
15X.12.lH/tI2.lH/. 12*//) 

X=1 

DO  100  I = I NDEX.LAST 
XT  AB ( F I  -SPACE  (  I  I 
100  F  =  F+1 

WRITE  (6.9CC1IN.X1.X2.X6ALL 

9  OOi  FORMAT (5X.17H  CABLE  DEFINED  BY.I5.15H  POINTS  FROM  X=.E15.8, 

1  6H  TO  X=,E15.8-20H  WHERE  X  AT  BALLOON1 .E 1 5. 8 ) 

F=1 

00  110  I = INDtX.LAST 
TTABIXI-STUFTI 1  ) 

110  F  =  F+1 

..'RITE  (6,9002)  DT.TTEST.DTPRT.DTMPRT.EPS 
9uC2  FORMAT  (19H  TIME  INCREMENT  DT  =  , E 12 . 5 ,4X , 1 3HT I  ME  TO  TEST- . E 12 . 5 .4X , 
115HDT  TC  PRINT  X 1  =  . E 1 2 . 5 . 4X , 24HDT  TO  PRINT  XI  AT  Ml  OPT - , 1 1 2 . 5  . /  .43 
2X.17H  TOLERANCE  ON  X I  - .  E15.8,/) 

F=  1 

DO  120  I = INDEX. LAST 
PTA6(X)=STUFP( I ) 

120  F  =  K  +  1 

WRITS  (6,9003)  G,D«E MO, DS.OuALL. OMEGA 
9u03  FORMAT ( 9X , 2f  G.18X.2H  D , 1 8X , 2HM0 , lbX , 2HDS . 13X , 12HS  OF  BALLOON, 11X, 
15HOMEGA,/,E17.8,5E20.8) 

C  SET  UP  CONSTANTS 

Cl  =  i:-T»DT/EMO 
C2  =  >  .  j«Cl/(DS*DS) 

C  3  *o , a*DT / EMO 

OXX-  .  5 v ( X  a  1 -XA2 I 

ALPHA  - 2 »  3025 85 1 / ( DXX+DXX ) 

J=3001 
NS=  1 

DO  15  0  1  =  1  ,N 
X=XTAB( I  I 
CAlL  RHOX(X.RHO) 

CALL  NU  (X.GNU) 

CALL  INTLRINVEL.NS»X.V*XV,VEL) 

R  =  A  I)  S  ( V  E  L  *  D*  P  T  At)  (  I  I /GNU  I 
CALL  CDR(R.CD) 

SPACE! J) ,0.5«RHO*D*CD»PTAB( I  I *PTAS( I  I/G 
PTAB!  I  I-VLL 
150  J=J+1 

C  VELOCITY  '.Ow  STORED  IN  PTAtS 

WRITE  1 1  ,  >  0  0  4 ) A  ,  X  A 1 , X  A 2 

9  004  FORMAT!  7X.4IH  V  TO  THc.  T = A . S I N ( OMEGA. T I . M X  )  WHtRt  A  =,tl5.8*l2H 
lBf.TwrEN  X  ",E15.8,8H  AND  X  =,E15.3) 

COUNT  = I  NO!  X .NMI D  -2 
SMID-SBAU  -DS+COUNT 
TMI')PT  =  T  I-Yn(NMID) 

C  PRINT  INPUT 

IMloEtM)  1 0 , 2  0 , 1 0 


}  /*■ 

? 

'  ^ 


•jf; 
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lu  Whi It  (6.900*1 

WR I T  t  I  6 . 9006 ) <  XV (1) . V ( 1 ) , 1  =  1 .NVELI 

9005  FORMAT(26X,28HlVELOCiTY  TABLE  ALT  I TUDE . 12X.8HVELOC1 T Y , / ) 

9006  FORMAT ( 37X.2E20. 3 > 

20  CONTINUE 

c  INITIALIZE 

T  =  -CT 
Txl=DTPRT 
T  TI D=DTMPRT 
TSTuP=2.0*TTEST 

c  initial  conditions 

DO  200  i*l.N 
X! ( 1.21=0.0 
200  X 1 ( 1.31=0.0 
X1M1DI 11*0.0 
M1DT0T=1 
MINI T=0 

C  THERE  hAvE  oEEN  M1DT0T  MIDPTS  STORED.  MINIT=0  MEANS  INITIAL  TIME 

OLDM1D=0.0 
IEND=0 
VALS  1  =  0 .0 
VALS2=0.0 

CALL  dOUNDIN.VAlS1.VAlS 2 1 

C.  oET  BOUNDARY  CONDI T IONS— X  I =0. AT  ENDS  OF  CABLE  SEGMENT 
C  RESET  LOOP 

3uo  T  =  T  +DT 

IF  (T-TSTOPI  301.800.800 
301  SVT  =A  •SIN(OMEGA*T) 

DO  250  1=1. N 
XI(I.1)=XI(I.2) 

2  5u  XI  ( I  ,2)=XI I  1 .3) 

C  COMPUTE  COEFFICIENTS  A.B.C.D  OF  XI  AT  N+l  SUCH  THAT 

C  A(M-l)  XMM-1J  ♦  BIM-ll  XI  I M)  +  C(M-l)  XKM+l)  =  D(M-l) 

CALL  AdCDIN. Cl. C2.C3. SVT. NXA1.NXA2.DT. ALPHA. XMID) 

C  GENERATE  ELEMENTS  IN  THE  SOLUTION  MATRIX 

CALL  GENERAT  INI 

C  SOLVE  THE  TRIANGULAR  SOLUTION  MATRIX 

CALL  SOLVE ( N 1 
ANSM 1 0  =  X  I (NMID.3) 

C  CHECK  THAT  DISPLACEMENT  IS  SMALL  NEAR  THE  ENDS 

NM  =  N-  1 

IF  I (AdSlXI (2*31 1+AbSCXl (NM.3) I 1-EPS)  350.350.340 
340  WRITE  (6,y010l  T.XII2.3) .XI (NM.3) .EPS 
9010  FORMAT  I 6H1AT  T=.E13.6.16H  XI  AT  THE  2ND  PT  =  . E 1 3 .6 . 22H  AND  XI  AT  TH 
IE  N-l  PT  =  »  E13.6.21H  WHtRE  THE  TOLERANCE* . E  1 3 • 6 ) 

WRITE  (6,9022)  T 

WHITE  (6.902  3! (XTAol I  )  .TTAol I )  ,XI I  I .31.1  =  1 »N  .NXPRT ) 

GO  TO  1 
35o  CONTINUE 

C  FIND  MAX, MIN  OF  TIME  HISTORY  OF  CABLE  MIDPOINT 

IF  (MINIT)  410,380,410 
380  MINI T=  1 

IF  ( ANSMIU-ULDMID)  400.390,390 
C  POSITIVE  SLOPE  FIND  MAX 

MPT=1 
uO  TO  500 

C  NEGATIVE  SLOPE  FIND  MIN 

400  MPT  =  - 1 

GO  TO  500 

410  If(MPT)  420,450.450 

C  MINIMUM  SEARCH 
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‘. c  „  ,f  I  ANSMIU -ULDMlN  )  300.AA0.4A0 
^  M )  N  )  MUM 

4  Uu  M  P  T  =  1 

WRITE  (6.9020)  ANSMI 0, T , SKI  0 . XMI 0 

9  _i2 1  FORMAT!/,  27H  MAXIMUM  DEFLECTION  AT  MIDPOINT - X I  =  ,E 15 . 6 ,4X , 

l&rtT IME=»E15.8,4X.2HS=El5.6.4X»2nX=,El5.6»/) 

GO  TO  500 

C  MAXIMUM  SEARCH 

ASu  I F ( ANSMI D-OL DM  ID)  A60.A60.50C 
C  MAXIMUM 

A60  MP1 =-l 

WRITE  (6,9021)ANSM1D.T.SMID,XMID 

902U  FORMAT!/, 37H  MINIMUM  DEFLECTION  AT  MIDPOINT - X i  = , E 15. 8, AX , 

15HTIME=,E1 5.3. AX. 2HS=E 15.6 ,AX.2HX=,E 15.8 ,/) 

50u  0LDM1N=ANSMID 

IF(T-TMID)  600.510,510 
510  TMlO=TMID+0TMPRT 
MI0T0T=M1DT0T*1 
XIMID(MIDT0T)=ANSM1D 
600  IF  ( ABS ( T-TX I )-0.01)  610. 610, 700 
610  T X I =  T X I +DTPRT 

WRITE  (6,9022)  T 

9 u22  FORMAT!//, 11X.6H  T I ME  =  E 15 . 8 , 7X , 2H  X.18X.2M  T.1BX.2HXU 
WRITE  (6,9023) (XTA6! I ) .TTAo! 1 ) ,X! ( I ,  3 )» 1  =  1 .N.NXPRT) 

9023  FORMAT(27X.3E20.8) 

700  IF  OEND)  770.710.750 
7 1 J  IF(T-TTEST)  300.720.720 
72u  MSUM=MPT 
1 END= 1 
GO  TO  3C0 

750  I F ( MPT  +MSUM )  300,760.300 
76 J  1 END=-1 

WRITE  (6,9022)  T 

WRITE  '6,90231 (XTAd( I ) .TTAo( 1  I ,X1  ( I »  3 ) » I  =  1 .N.NXPRT ) 

MSLM  MPT 

Hj  » F  (MPT+MSUM)  300.800,300 

3uw  WR )  )  E  (6,9022)  T  , 

WRITE  (6,9023) (XTAB( 1 ) .TTAbl I ) ,X1 ( 1 ,3) ») =1  .N.NXPRT) 

WRITE  (6,9025)  SMID.XMID.TMIOPT 

/ 2  3  FORMAT  (  1HI ,  13X.39iiTIME  itlSTCRY  OF  CABLE  MIDPOINT  WMcRt  S«.tla.o, 
lJX<.2rlX  =  ,Eib.b«iX,2hT-»Eji3»b</  ,49X,2bT  ,  lbX,2H>  )  ,/) 

T 1  ME  =0.0 

00  900  1 = 1 ,M 1 DTOT 

WRITE  (6,5026)  T1ME.XIM10I I  ) 

■  .26  FORMAT (37X.2E20. 8) 

90u  TIME=TIME*DTMPRT 
WRITE  (6,9030) 

;J3u  FORMAT (// .A9X.22H  #****END  OF  CASE*****) 

GO  TO  1 
99  CALL  EXIT 
END 

SUBROUTINE  ABCDIN.  Cl ,C2 .C 3 .SVT »NXA 1 , NX A2 ,DT .ALPHA . XM) D ) 
DIMENSION  SPACE! 100001.XI ( 2000 . 3 ) , T T Ab i 2000 ) .PTAl3(2000) 

1 ,X  T  AB ( 2000 ) 

COMMON  SPACE. XI  .TTAB.PTAb.  IOtcM.  )OEbd»  iL’EeA,  lutbG,  )OLbS 
1 . X  T  Ab 

c  computes  n-2  colficients  a.o.c.ano  o  ott-iNiNu  xi  at  n+i 

c  olVLN  XI  AT  N  AND  AT  N-l 

C  SPACE! 1-2000)=A,  ( 200 1 -4COO  )  =b ,  ( 4CC1-6000  I  =C.  ( 6  00 1 -cooo ) =D 

XI (M,  1  ) =  X I  (M.N-1) .XI (M.2 )  =  x I (M.N) 

HM  -N-  1 


!>»£*:»  .-T- 


•# 


jef 


4 

*M. 

t: 


DO  1  CO  M=2,NM 

MM=M-1 

MP=M+1 

SlOP=SPACE  (OCC0*M) 

C  ADJUST  FORCING  FUNCTION 

VEL=PTAB(M) 

VV  =  VEi_*ABS(  VEL) 

IF  (M-NXAl)  6C.50.40 
40  I F  (  M-NXA2)  50*50*60 
50  Dxx=XTAd(M)-XM10 

POWER  =  -AlPha*DXX*DXX 
C  SVT  = A. S I N IOMEGA*  T ) 

VT-SvT»tXP(POWER) 

C0N=VEL+VT 
GO  TO  70 
60  C0N= v£L 

7o  CONAys=AHS(CON-(Xl (M.2I-XI (M.l) »/OT) 

CONST -C3‘SLCP*CONABS 
A=-ITAo(MK)*C2 

0=1 .0*1 1 TAb(M)+TTAb(MM) >*C2+CONST 
C  =  - T  T  A3 ( M ) *C2 

D  =  -Xl(M.l)-*2.*XI(K.2)-C»(2.0*XI (MP,2)+XI IMP , 1 ) -2 . 0*X I !M,2l 
1-Xl (M.l) )+A*(2.0*XI (M,2)+XI (M.l)-2.0*Xi(MM,2)-Xl (MM»1)  ) 

2  +CONST»Xl (M.l)  -C1«SL0P*(VV-C0N«C0NA6S) 

SPAC£(MM)=A 
SPACE  (  MV, +  2000  )  =  B 
SPACE (MM+4C00)-C 
l(ju  SP*Ct (MM+6000) =D 
RETURN 
END 

Subroutine  generat  in> 

DIMENSION  SPACE (  10000) *x 1(2000.3) .1  TAB ( 2000 ) *PTAB( 2000 ) 

1  *  X  T  A3 ( 2000 ) 

COMMON  SPACE. Xl *TTAd.PTAo*IDEbM.IDEbB.IDtbA,IDEBG.IDtBS 
1.XTA3 

u  ucNcRATr.  r.LcMcNTS  IN  Tfit  SOLUTION  MATRIX 

C  SPACE! 1-200C)=A,  (2001-40001=3.  ( 4001-6000 ) =C.  ( 600 1-6000 ) =-D 

C  GtNERATtS  C(  D-CIN-l)  ,DU)-D(N-1  )  ,, 

SPACE ( 4  „C  1 )= SPACE! 4001) /SPACE (2001) 

SPACE (6w01)=SPACE(6C01)/SPACE(2001) 

NV=N-3 

DO  I  00  M  = 1  * NM 
MP=M+1 

M2000  =  MP  »  2000 
M4 JOc=MP*4000 
M  6  0  0  c  =  MP  +  6  0  0  0 

SPACElM2C.^O)=SPACt(M2000)-SPACt(MP)XSPACE(M+4000) 

SPACE (V4G 00 ) =SPACE (M400u) /SPACE (M2C00  > 

IuJ  SPACE ( M60CO )  =  ( SPACE ( M6000 ) -SPACE (MP ) *SPACE ( M  +  6000 ) ) /SPACE  (  M2000 5 
RETURN 
END 

oUJRct'TlNE  SOLVE  (N) 

J I  Me.  NS  I  on  SPACE!  1 1 000  )  .XI  (2000. 3  >  .TTABI  2000)  .PTAo(2000) 

1  ,XT  Af3(  2000) 

COMMON  5  PACE . X  I . T T AB . P T A d , I DtoM . IDEoB. I  DcbA  ,  1 DEBG *  1 DlBS 
1  .  X  T  A  B 

c  SOLVE  T i  E  TRIANGULAR  SOLUTION  MATRIX 

c  SPAClI 1-2^00  )=A,  (2001-40C0) =B.  ( 400 1 -60 00 ) =C  .  ( 600  1-O0G0 ) =D 

NM-N-2 

1 NDcX  =  NM  +  600  0 

XI  IN-  l»3)=3PACt(  INuCX) 
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00  100  1=2. NM 
J=NM  -1+1 

100  XI ( J.1.3)=SPAC£< J  +  6000) -SPACE <  J+40C0t«Xl ( J+2,3) 

RETURN 

END 

subroutine  bound  <n.vAi_si.vals2) 

DIMENSION  SPACE! 10000 ) .XI ( 2000. 3 1 . T T Ad ( 2000) .PTABI2000) 
1.XTABI2000) 

COMMON  SPACE. XI , T TAB .PTAB . 1 DEbM. 1 DEbB. I DEbA, 1 DtbG » I  DEBS 
1 .XT  AB 

C  SET  BOUNDARY  CONDITIONS — Xl =0 »DX I /DT*0  AT  ENDS  OF  CABLE  SEGMENT 
Xl ( 1.3I=VALS1 
Xl ( N. 3 ) =VAL52 
RETURN 
END 

SUBROUTINE  NU  (X.-GNUl 
C  K.INEMAT  I C  VISCOSITY  VS  ALTITUDE 
I.  (X-10000.)  1.1.2 

1  GnU=-4.823+X/26800. 

GO  TO  5 

2  IF  (X-17700.I  3.3.4 

3  GN'J=X/  17100. -5.035 
GO  TO  5 

4  GNU-X/ 14250. -5.2421 

5  GNU=10.**GNU 
RETURN 

END 

SUBROUTINE  INTER  i N.NS.Xl .V ,XV ,SVL ) 

C  GIVEN  THE  ALTITUDE.  THIS  ROUTINE  INTERPOLATES  LlNtARLY  TO 

C  OBTAIN  THE  *1ND  VELOCITY  FROM  INPUT  TABLE 
DIMENSION  V(  1001  .XV! 100) 

NN=N-NS+1 

IF  (XL-XV(NN))  5.5,6 

6  NS  =  1 

5  DO  1  I=NS.N 
J=N-I+1 

IF  (XL-XV(U))  1.3.4 

3  SVL=V( J) 

GO  TO  2 

4  SVL  =  ( V ( J+ 1 ) -V ( J) ) *(XL-XV( J+l )  )/ (XV! J+l )-XV( J ) )+V( J+l ) 

GO  TO  2 

1  CONTINUE 

2  N  S  =  J  +  1 
RETURN 
END 

SUBROUTINE  CDR(R.CD) 

C  DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER  FOR  CABLE 

IF  (R-2.23)  1.1,2 

1  CD=10.8*R**<-.742) 

GO  TO  10 

2  IF  (R-8.0)  3.3,4 

3  CD=9.15*R**<-.526) 

GO  TO  10 

4  IF  (R-1000.0)  5.5,6 

5  CD=4.95*R**(-.232) 

GO  TO  10 

6  IF  (R-lOOOO.)  7,7.8 

7  CD= 1.0 
GO  TO  10 

8  CD=1 .15 
Is  RETURN 


71 


,  Wrn&imammxum*.*  *> 


END 

SUBROUTINE  RHOX(X.RHO) 

wass  density  vs  altitude 

RhO=.25-X/160DO. 

SHO=lC.»««rtO 

RETURN 

END 


t 
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APPENDIX  D 

WIND  INDUCED  VIBRATIONS  BY  MEANS  OF 
VORTEX  SHEDDINGS 


V 


Considering  the  cable  subjected  to  a  steady  wind  flow,  any  high  frequency 
dynamic  effect  on  the  cable  will  cause  it  to  vibrate  about  its  static  con¬ 
figuration.  In  other  words,  the  dynamic  response  of  the  cable  to  high 
frequency  excitations  can  be  considered  as  a  perturbation  on  a  static 
initial  configuration  which  is  related  to  a  steady  wind.  Further  simpli¬ 
fication  of  the  problem  can  be  done  by  considering  the  high  tension  in 
the  cable  and  its  low  weight  per  unit  length.  The  propagatv^,;  velocity 
of  the  transverse  wave  in  the  cable  is 


_ _ 2 


c  ■  \/S 


In  addition  the  frequency  of  shedding  vortex  pairs  is 

sT(R)voP 

v  '  D 


(D-2) 


where  S^,  is  the  well  known  Strouhal  number.  In  the  present  case  f^ 
is  of  the  order  of  1000  cps.  The  wavelength  which  one  should  be  con¬ 
cerned  with  is 


2L 


C_ 

f 

v 


f 

4 


t 


which,  for  the  present  purposes,  is  of  the  order  of  2  meters.  '  once 
the  most  effective  mode  of  the  cable  will  be  the  one  which  is  of  wave¬ 
length  2L. 

In  order  to  estimate  the  response  of  the  cable  to  vortex  sheddings,  one 
can  solve  the  equation  of  motion  for  the  cable  considering  a  portion  which 
has  the  length  L.  This  portion  can  be  assumed  to  be  supported  at  its 
e ndpoint  s . 
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The  perturbed  dynamic  equation  of  motion  for  the  cable  is 


-  T  pCx.VlPZD  sin2rrf  t 
2  H  K  o  v 


(D-  3) 


where  is  the  lift  coefficient  (Refs.  5  and  6)  and  V^P  is  the  normal 

component  of  the  wind  velocity  vector  with  respect  to  the  cable. 


The  last  equation,  which  has  been  solved  numerically  by  means  of  finite 
differences  is  discussed  in  a  following  section.  The  scheme  that  was 
used  for  this  purpose  was  an  unconditionally  stable  one.  The  integration 
process  ran  on  a  computer  for  two  hundred  cycles  of  the  forcing  function. 
Within  this  process  of  integration  the  cable  reached  a  steady  state  re¬ 
sponse.  For  L  =  0.775  meters  the  maximum  deflection  of  the  cable  was 
_4 

0.58  x  10  meter,  which  does  not  indicate  any  considerable  increase 
in  the  stresses  in  the  cable.*  In  conclusion,  the  cable  subjected  to  high 
tension  converts,  by  means  of  the  vortex  sheddings,  wind  flow  energy 
to  other  kinds  of  energy;  however,  there  is  no  indication  that  the  high 
frequency  vibrations  of  the  cable  with  small  amplitude  might  cause  any 
failure. 


*The  properties  of  the  cable  and  wind  field  were  chosen  from  the 
steady-state  example  and  are  as  follows;  midpoint  of  cable  portion  at 
4000  m;  T  =  773.0  Kgf,  P  =  0.522,  P  -  1.0  Kg/m3,  y0  =  13.8  m/sec, 
m  =  8  x  10-4  Kg/m,  D  =  2.  5  x  10~3  nr,  Cq  =  1,  =  0.  76,  =  0.  22. 
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Numerical  Analysis  of  Problem 


Equation  of  Motion 


_b  ( rp  iW  N  _ 

3S  \  m  .  -2/  mo 


™  ~  DV2  P2  V*  C 

^2  2g  o  m  K 


+  TT  CD  (Rm)  ° 


P2  + 
m 


/  \2~11/2 
(5W\  d_W 

\dt  /  bt 


where  W  is  a  deflection  i  to  the  plane  of  Fig.  1.  The 
short  cable  section  length  leads  to  assuming  the  following 
parameters  constant:  V T m>  Pm>  p(xm)  ,  CD  (Rj  , 

CK  (Rm)  '  ST  (Rm) 

The  equation  of  motion  becomes: 


-  -£-C„D  V2P2Vt  + 


o  at2  2g  K  o 

•[V  ♦  (f)T 


2g  D 


whe  re : 


V1  =  sm27Tfvt 


S,r  V  P 
1  o 


given  the  frequency  f^,  the  following  parameters  can  be 
calculated: 


At 


?0f 


-t  = 


2f 


max 


j—*  *(no.  cycles) 
v 


b)  Difference  Analogue 
T 


4  AS 
+  W 


(" 


,  1W  XI  x,  +  W  ,  X,  ■  2W  ,,  +  2 

2  1  m+l,n+l  m-l,n+l  m,n+l 


.  -  2W 

m-  1 ,  n  m 


.n) 


+  W  .  +  W  .  .  -  2W 

m+ 1 ,  n-  1  m- 1 ,  n-  1 


(w  4.1 

\  m+l,n 

^m,  n-  lj 


m  / 

=  _o  fW 

At2  V  m.n+l 


2W  +  W 
m,  n  m 


\  -  P 

,  n-lj  2g 


-  P  Cr.  D  V2  P2  V1 


pcdd 

2g 


2  2 
V  P4  + 
o 


(w  -  W 
'  m,  n _ m,  n-1  / 


pM? 


At 


C 


w  ,  -  w 

m,  n+1 _ mj 


2  At 


*D 


The  cable  segment  is  analyzed  for  half  of  its  length  t  and 
is  defined  between  points  1  and  M  where  M  is  at  1 /2, 
The  "ficticious"  point  M  +  1  is  used. 


c) 


Initial  Conditions 


W  -  0.0  for  m  -  1,2,  .  . .  f  Ivi  +  1 

m,  o 


d)  Method  of  Solution 


If  A 

in-  l 

W  4 

m  -  1 ,  n+ 1 

Bm-1 

=  D 

m-  1 

for  rn  -  2.  . 

.  .,  M 

+  C 


W 


then: 


At"T 

4AS2m 


At  2  2 

=  1  -  2A  .  +  — - —  VZ  P** 

m- 1  4m  g  o 


(w  -  W 

\  m,  n  m,  n-  1/ 


=  -  A  ,  W  A1  .  +  W  .  -2W 

m-1  m+1,  n-1  m-  1 ,  n- 1  m,  n- 1 


t  fw 

m-  1  I  m 

(W  X, 

'  m+l,n 


+  W  -  2W 

m-1,  n  m 


..)] 


At2  pC„  D  V2  P2  V1  At  pC  D 

+  K  °  +  D  yZPZ 

2m  g  4m  g  o 


fw 

\  m,  n 


yT 

i,  n-  1/ 


m,  n- 1 


+  2W  -  W 

m,  n  m ,  n - 1 


The  solution  of  the  linear  system  lor  W  is  carried  out 

'  m,  n 

by  the  matrix  inversion  procedure  described  in  Appendix  C. 
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e)  Boundary  Conditions 


Wl,n+1 

=  0.0 

AM-1 

-  am-i  + 

CM-1 

dm-i 

=  DM.l  + 

CM..'  [WJ 

+  Z(WM+l,n  -  W 

-  W, 


M- 1 ,  n- 1 


f)  Solution  of  System  of  Equations 

If  the  following  computations  are  made  initially. 
Cl  =  Cl/Bl  D1  =  Dl/Bl 


The  generating  sequences  maybe  given  for  m  =  1 . 


M  -  1 

3m  +  l 

- 

Bm+1 

'  Am+lCm 

Cm+1 

= 

Cm+1 

Bm+1 

Dm  +  1 

Dm+1 

"  Am+lDm 

Bm  +  1 

F  ollow 

ing 

the  computations  of  the  above  sequences  the 

solution  for  W(m,t)  may  be  written  for  m  =  M+ 1 ,  M 

Z 

W 

M,n+1 

- 

dm-i 

W 

M  +  l , 

n+1 

= 

WM- 1 ,  n  +  l  '  2(WM+!,n  ‘ 

_  W  +  W 

M+l.n-l  M-  1 ,  n- 1 


«y  "1 ' 


-iPiaww-w  -1  »,j“ 


j 


w 


=  D  ,  -  c  ,  W  _L1 

m-1  m-1  m+l,n+l 


m,  n+ 1 

where  m  =  M-i,  M-2,  ...,  2 


g)  Required  Data 


Input  constants:  D,  mQ,  g,  Vq,  T,  P,  p,  CD>  C  ST 
whe  re  CD’  CK*  and  are  functions  of  Reynolds  No.  R 


I 
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W 


Input  to  Program  MUSIC 


Columns  10 

20 

30 

40 

50 

60 

Card  1 

NDA 

NMO 

NYR 

■ 

Card  2 

NPOINT 

NCYCLE 

NXPRT 

NTPRT 

IDEBUG 

B 

Card  3 

G 

X 

DS 

T 

P 

RHO 

Card  4 

EMO 

D 

CD 

CK 

ST 

Notation 

NDA,  NMO, 
NYR 

NPOINT 


NCYCLE 

NXPRT 


NTPRT 

IDE BUG 


G 

X 

DS 

T 


Date 

Number  of  points  in  the  cable  section  to  be  analyzed. 
The  section  will  be  analyzed  up  to  its  midpoint 
M  =  (NPOINT /2)  +  1 

Number  of  cycles  to  run  program 

Increment  used  in  printing  deflections  W  at  any  time 
(i.e.,  printed  at  every  NXPRT  point  out  of  M  +  1 
points) 

Increment  used  in  printing  the  time  history  of  the  mid¬ 
point  (i.e.,  printed  at  every  NTPRT  time) 

If  IDEBUG  =  1  deflections  at  every  NTPRT  time  will 
be  printed 

If  IDEBUG  =  0  only  the  time  history  of  the  midpoint 
will  be  printed 

Gravity  constant  (meters/secL) 

Altitude  oi  the  cable  section  (meters) 

AS  =  S  increment  between  the  points  in  the  cable 
section  (meters) 

Tension  in  the  cable  section  (  kilograms-weight) 


m: 


i 

{ 
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P  P  =  Ax/  AS  =  cos  a  =  determines  slope  of  the  cable 
section 

RHO  O  =  density  of  atmosphere  (the  program  divides  p  by  g) 
(Kg/meters^) 

VO  =  wind  velocity  on  cable  (meters/sec) 

EMO  M  =  mass/unit  length  of  cable  (Kg/rneter) 
o  n 

D  Diameter  of  the  cable  (meters) 

*\ 

Op  Drag  coefficient 

Lift  coefficient  \  functions  of  Reynolds  number 

S,j,  Strouhal  number 
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* 


,_v: 


r\  r»  r. 


upoGRAM  music 

dynamic  ANALYSIS  OF  THt  MUSICAL  PROBLEM  SN3AA-716 
DI MtNSI ON  *(100.3).  A(  lliu ) » B  < 100) » C  I 100) » D( 100) 

COMMON  W,A,=,c,0 

NPOINT  POINTS  OtSCRlBt  A  SHORT  StGMfcNT  OF  THt  BALLOON  CAdLE 
C  THt  SEGMENT  IS  ANALYZED  FOR  M=NP01NT/2+l  POINTS  USING  A  FICTITIOUS 
C  POINT  1 

C  THt  ShORT  StGMfcNT  ASSUMtS  CONSTANT  VALUES  OF  X»TiP,kHO*VU»CD»ClC»ST 

1  READ  (S.eLGUNDA.NMO.NYR 
8  0uu  FORMAT (  7110) 

I r  (NOA)  99.99*2 

2  READ  (5.8CC0)  NPOINT .NCYCLE .NXPK T *N TPRT  .  I  DEBUG 

READ  IS. 8001)  G.X.OS.  TtN  .P  .RHO  •  VO.  EMO*  01  A. CO.  CIC*ST 

cjOI  FORMAT  (7F1C.8) 

C  TMAX  OF  THE  RUN  IS  CALCuLATtO  AS  NC YCLt . 2P I /FKtOUtNCY  FV 

C  DT  OF  TnE  RuN  IS  CALCuLA TtD  AS  1/(20. FV) 

■  HE  CABLE  SEGMENT  LENGTH  1b  CALCULATED  AS  SORT  I T/MO ) / C 2. FV ) 

DATA  IS  PRINTED  AT  tVERY  NXPRT  POINT  FOR  tVERY  DTPRT  TIME 
DTPRT=2PI/(FV.NTPRT) 

WRITE  16,9000)  NDA.NMO.NYR 
FORMAT  1 1 Hi ,26X, SOHO YN AM  I  C  ANALYSIS  OF  THt  MUSICAL  PROBLEM— SN34M 
1716  ,4X , I  2 , lH/» 12. 1H/ .12,//) 

Compute  constants 
cy.:lls=ncycl.e 

F/=ABS<VG*ST*P/0I A) 

0 T-..0S/FV 
TP=NTPRT 

DT^(a-6.2631653/(FV*TP) 

CA=LK=  C.S«SORT(TEN/EM0)/FV 

OMEGA-6. ;oJlb53«FV 

TMAX  =  6.2o318S3*CYCLES/FV 

DT2=DT»CT 

DS2=DS»DS 

Cl=  v.2S«DT2*TEN/(tMO*DS2) 

C2=  0.2SvRf'C*CD*DIA*DT/(EM0*G) 

C3=(VO'iP)**2 

PHI  Is  r;.S*RnO«C<*DI  A*C3/G 
C4  r  Ph  I 1*DT2/EM0 
POINT/2  +1 
.'iP  ;  M  +  1 

I  It  16.9001)  NPOINT, M. CABLE. X »RHO . VO .EM0.D1 A. TEN.P 
9^1  F'uRM.ATl  19X.25H  CABLE  SEGMENT  DEFINED  BY.13.24H  POINTS  AND  ANALYZED 
1  FDR.I3.26H  POINTS  UP  TO  THE  MI DPOI NT . // » 23H  CABLt  DATA  LtNGTH 
2-.E15.S.5H  X=»E15.8,1SH  DENSITY  RHO=  ,t 15 , 8 » 15H  VELOCITY  V0=. 

3LlS.b»/,lSX»9H  MASS  M0=»E1S.8.  13H  DIAMETER  0-.E15.6,  12h  TtNSIO 
<•:,  Tr.ClS.8,  10H  PrDX/DS  =  »ElS«8./) 

...RITE  I  6  ,  V  0  0  2  )  CD.CX.ST 

,uu2  f orm a t i i 9h  functions  of  Kt ynolds  number, isx,3hc  =,tis.», lsx, jhc  =  , 

ItU'.o, 1SX.3HS  -,E 1S.8./.AAX.2H  D » 3 2 X , 1HK , 3 2X , 1 HT , / / ) 
wKI TL I 6,9«0i  )  TMAX.OT .DTPRT. OS.G.FV 

/■-.03  FORMA i ( 3X. 13H  MAXIMUM  T I  ME . 8 X  .  1  3H T I  ME  1 NTtRVAL » 4X . I 7hPR I  NT  I NG  I NTE 
1RVAL .8X.7HCELTA  5  .  13X,7hGRAVI TY » 1 IX  .  12ilFRtOUENCY  FV»/»tl7.8,5E20.H 
2) 

initialize 

7  -  -DT 

TPkINT-DTPKT 
Du  UO  1  =  1. MP 
. I  1 .21=0.0 
iL^.  X  (  I  »  I  )  =  u  •  0 

..RITE  (  6  ,  *0 1 5  ) 

» U  I  3  I  U  •.’■’.AT  I  1 1! I  .  AOX  «  30HT  IME  HISTORY  OF  LAbLt  SEGMtNT  M 1  DPO I N T  » / /  » *■  * k  * 


15m  TlME.iax. 1HW.13X. 13H2ND  QER  Dw/DS • 1 IX * 5MDW/DT > 

LOCPs-l 
U  1  S  T  =  1 

C  RESET  LOOP 

2 Ov  LOOP -L JCP+ 1 
SLGGP=LOCP 
t=slcop*ot 

VT=SIN(0M£GA*T 1 
DO  300  1  *  1 »MP 
*11  •  1 )  - * (  1«2) 

30u  w( I ,2 >  =W( I .3  ) 

C  COMPUTE  COEFFICIENTS  A.b.C.D  OF  w  SUCH  THAT 

C  A(M-l)  W(M-l)  +  BlM-1)  W(V.j  ♦C(M-l)  W(M+1  )=D(M-1  ) 

CALL  A3CD(M,C1.C2.C3.C4,VT.DT) 

C  ADJUST  BOUNDARY  CONDITIONS 

CALL  EOUNO  (Ml 

C  GENERATE  ELEMENTS  IN  THE  SOLUTION  MATRIX 

CALL  GENERAT(M) 

C  SOLVE  THE  SOLUTION  MATRIX 

CALL  SOLVE  i  K.D2WDS2  >  DwDT  «  Db2  «DT  ) 

IF  (IDESUG)  45C.45C. 350 
350  IF  (T-TPRINTI450.400.400 
40l  PMI=PhI1*VT 

WRITE  (6,9010)  T.PHI 

9010  FORMAT!  87X,  6H  T I  ME - • 1 10. 3 . 3X .4HPH I  =  , 1 10. 3 . / . 97X ,3H  PT.13X.1HW) 
WRITE  (6,9011) (!<*(!. 3)*lsi.M«NXPRT) 

9011  FORMAT  ( 96X . 1 4 , E2 0 .8  ) 

450  IF  (T-TPRINT)  500.460.460 
46 J  L IST=LI ST+1 
SL I S  T  =  L I  ST 
TPRINT-SlIST»DTPRT 

WRITE  (6,9016)  T.W5M.3) .D2wDS2.DwDT 

9C  16  FORMAT! 17X.4E20.8) 

500  IF  (T-TKAX)  200.600.600 
60 J  WRITE  (6,9020) 

9u2w  FORMAT!//. 49X.22H  *»***EhlD  OF  CASE*****) 

GO  TO  1 

99  call  exit 

END 

SUbROU T I Nt  ABCD(M,C1.C2.C3.C4.VT.DT) 

DIMENSION  '*(100,3),  A ( luo )  , B  1  100 )  , C <  1 00  )  , D (  100 ) 

COMMON  w, A.B.C.D 
DO  ICO  1=2, M 
MM- I -1 
MP  =1+1 
A ( MM ) =  -C 1 

CON  =  C2»SCRT IC3+ ( ( (W< I ,2)-W(  1 ,1 1  ) /DT  >**2)  ) 

B ( MM ) =  1 . C- 2 . 0*  A ( MM | +  CON 
C (MM) =A(MM) 

DIMM) =-A(MM) » I  WIMP, 1 ) +W (MM, 1 )-2.0*W( I , 1 ) +2.0*1 W(MP,2 ) +w(MM,2) 

1 -2.0*W( I ,2) ) I ♦L4*vt+CON*W( I ,1  )  +  2.0*  '  ( 1,2) — w (1.1) 

lOv  continue 

RETURN 

lND 

subroutine  BOUND ( M ) 

DIMENSION  W (  100,3)  ,  A( 10o)  ,B( 100) ,C( 100) ,D( 1001 

COMMON  W.A.B.C.D 

•'11,31=0.0 

MM=M-1 

M  P  =  M  ♦  1 

A  I  MM )  =  A  (  VM  )  +  C  (  MM ) 
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L 


* 


0 ( KM  1  -C (KKI +CIMM) * (WIMP* 1 ) -W(MM. 1 ) +2.0* ( W(MP*2 ) “W (MM *2 ) )  ) 

Return 

eno 

subroutine  generatiM) 

DIMENSION  W ( 100 1 3 1 •  A(  100) .  B  ( 100) »C  ( 100) *  D( 100) 

Common  w.a.b.c.d 

C<  1  )  =C  <  1  )/B(l) 

D(1)*B<1)/8<1) 

MH2*M-2 

DO  ICO  1=1. MM2 
MP=I+i 

b  ( KP  )  =  B  ( XP  )  -  A  (  MP  )  *C  (  1  ) 

C  ( MP  )  =  C  <  MP  )  /  B  ( MP  ) 

D(KP)=(D(MP)-A(MP)*D(I n/b(MP» 

10U  continue 

RETURN 

END 

SUBROUTINE  SOLVE (M,02W0S2.DwDT.DS2. DTI 
DIMENSION  W ( 100.3)  •  A(IOO). 8(100). C(IOO). 0(100) 

COMMON  A.AtB.C.D 

MM  =  M-1 

MP=M+1 

W(M.3)=0(MK) 

DO  100  1=2. MM 
J=MM- 1*2 

*(J.3)=D(J-1)-C(J-1)*W(J+1,3) 

100  CONTINUE 

DwOT  =  0.5*(VMM.3)-w!K,l))/DT 

SYMMETRICAL  CONDITIONS  AT  POINTS  M-l  AND  M+l 

W I  MP .  3  )  =  ft1  ( MM .  3 )  -2 . 0*  ( W  ( MP .  2 )  -W  ( MM.  .2 )  )  -*  l  MP  » 1 )  +*(  I  MM  .  1 ) 

D2*D£2  =  C.25*(W(MP.3)+W(MM,3)  +  W<MP,1 ) +W < MM . 1 ) +2 . 0* ( W ( MP  ,2 ) +W ( Mrt,  2 ) 
1  -w(M.3)-W(M, 11-2.0*  W1M.2 ) ) I/DS2 
RETURN 
END 
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